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^methods  can  also  be  used  to  solve  nonasymptotically  stable  differential  equations  since 
they  are  a  generalization  of  Linear  Multistep  (LMS)  methods. 

..  i 

This  report  presents  a  detailed  formulation  of  NLMS  methods  and  discusses  a 
FORTRAN  V  computer  package  specifically  developed  to  implement  NLMS  methods. 

The  package,  also  available  in  ANSI  FORTRAN,  is  presently  operational  on  Univac  1108, 
IBM  360/370,  and  CDC  6600  computers.  Several  desirable  features  are  included  in  the 
computer  program,  namely,  a  fixed  or  variable  step  size,  self-start,  a  selection  of 
characteristic  polynominal  coefficients,  predict- and-correct  m  times,  and 
the  inclusion  of  LMS  methods.  ^Whenever  matrix  A  is  a  function  of  t,  a  periodic  de¬ 
composition  technique  is  emplltecl.  Otherwise,  a  decomposition  technique  can  be 
selected  by  the  user.  Nonlineip  mqltistep  methods  and  their  associated  features 
constitute  a  powerful  means  for  solwng  initial  value  problems  of  stiff,  nonstiff,  linear, 
and  nonlinear  ordinary  differential  equations.  The  validity  of  NLMS  methods  has  been 
proved  theoretically  while  the  effectiveness  of  the  computer  program  will  be  supported 
by  the  numerical  evidence  to  be  presented. 
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NUMERICAL  SOLUTIONS  OF  INITIAL  VALUE  PROBLEMS 


1.  INTRODUCTION 

To  solve  nonstiff  differential  equations,  the  conventional  Runge-Kutta 
and  linear  multistep  methods  are  typically  employed.  These  conventional 
methods  are  impractical  for  solving  stiff  equations  because  prohibitively 
small  step  sizes  are  required  for  accuracy.  To  overcome  this  difficulty, 
nonlinear  multistep  (NLMS)  methods  can  be  applied;  NLMS  methods  have 
been  shown  to  have  advantages  over  existing  stiff  methods.  It  is  desirable 
to  be  prepared  with  a  variety  of  effective  methods  for  handling  both  stiff 
and  nonstiff  equations.  The  search  for  such  an  efficient  program  will 
probably  never  end.  To  be  effective,  a  package  must  be  reliable,  easy 
to  modify,  and  convenient  to  use.  In  addition,  other  features  are 
desirable.  To  satisfy  these  needs,  the  NLMS  method  was  developed.  It  is 
now  known  that  strongly  stable  NLMS  methods  are  consistent  and,  therefore, 
convergent.  The  complete  theory  of  NLMS  methods  has  been  published. 1 

This  report  will  first  outline  some  preliminary  considerations  and 
assumptions  and  then  describe  the  formulation  of  NLMS  methods  along  with 
the  various  built-in  features.  Since  the  principal  objective  was  to  solve 
stiff  equations,  variable-order  techniques  were  not  utilized  because  NLMS 
methods  of  all  orders  were  considered  to  be  effective. 

The  various  features  are  described  and  illustrated  by  problems. 

A  section  of  numerical  results  describes  the  selection  of  these  problems  from 
well-known  sources,  categorizes  them  In  various  classes,  and  then 
summarizes  their  characteristics  in  tabular  form.  Various  NLMS  methods 
are  applied  to  solve  these  problems  with  given  initial  values.  The  computed 
solutions  are  next  compared  with  exact  solutions.  The  computational  accuracy 
is  measured  mainly  by  relative  error,  which  will  be  defined  in  the  next  sec¬ 
tion.  An  error  of  10-*0  is  used  as  an  upper  bound  to  acceptable  performance 
in  test  examples.  Numerical  results  are  printed  out  to  eight  significant 
digits.  The  User's  Guide  is  intended  to  help  the  user  become  familiar  with 
the  Inputs.  To  assist  in  this  objective,  the  way  to  start  a  problem  with 
various  sample  setups  is  described  within  each  test  problem.  Only 
FORTRAN  V  inputs  are  described  since  ANSI  FORTRAN  Inputs  are  made 
optional  to  the  user.  The  usage  of  a  few  user-supplied  or  optional  sub¬ 
routines  is  also  described  within  each  test  problem,  which  exercises  various 


1 


TR  5341 


* 


options.  Programs  were  originally  written  in  FORTRAN  V  language  and 
checked  out  on  the  Univac  1108  computer  using  double-precision  arithmetic 
by  means  of  the  EXEC  8  operating  system. 

These  programs  are  also  available  on  IBM  360/370  and  CDC  6600  in 
FORTRAN.  A  version  of  the  program  is  made  available  In  ANSI  FORTRAN, 
where  the  external  and  the  adjustable  dimensions  were  purposely  eliminated. 
Detailed  operational  descriptions  of  the  IBM  and  CDC  machines  are  not  in¬ 
cluded.  In  the  appendix,  a  program  listing  of  both  FORTRAN  V  and  ANSI 
FORTRAN  are  given. 
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2.  PRELIMINARY  CONSIDERATIONS 

In  this  section,  we  define  the  problems  under  consideration,  give  the 
norm  used  for  comparison,  and  state  some  assumptions. 

1.  Let  us  consider  the  initial  value  problem  of  a  system  of  first- 
order  ordinary  differential  equations  of  the  form 

y' =f(t,y);  y(tQ)=y0. 

which  can  also  be  written  as 

y'  “  Ay  +  g(t,y)  ;  y<V=y0 

in  the  region  R,  defined  by  -«°<a  <t<b<«  ;  ||y||<«o.  The  matrix  A  is 

nonsingular,  either  a  constant  or  a  function  of  t.  We  assume  that  our  initial 
value  problems  satisfy  the  conditions  required  by  the  existence  and  uniqueness 

theorem;  therefore,  f  and  9  both  satisfy  a  Lipschitz  condition  with  Lipschitz 

constants  L*  and  L,  respectively,  and  f  is  continuous  in  R. 

2.  Subject  to  the  consideration  1.,  we  consider  only  first-order 
equations.  This  is  not  a  restriction  because  higher  order  equations  can  always 
be  decomposed  into  a  system  of  first-order  equations.  Equations  that  are 
stiff,  nonstiff,  linear,  nonlinear,  homogeneous,  and  nonhomogeneous  are  all 
taken  into  consideration.  We  regard  every  problem  as  a  system  of  equations. 

3.  Let  yc  be  a  computed  solution  vector  andy£  be  an  exact  solution 
vector.  The  norm  used  here  to  measure  the  error  is  defined  as  follows: 

Hyli’p^Jlyllp  =  T  lyjl*  p-3) 


(2.1) 


(2.2) 


We  define  the  relative  error  E  to  mean 


E  = 


max 

i 


(yt>  -  (y4) 

c _ _ 


(2.4) 
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3.  FORMULATION 


The  linear  multistep  methods  of  step  K  can  be  expressed  by 
K  K 


E  rn+i =hE/,  *n+i  • 

i=0  i=0 


(3.1) 


where  +0,  I  <*  I  +  |  fl  I  >  0,  and  and  8,  are  constants  Independent  of  mesh 
K  o  o  i  ri 

size  h.  An  LMS  method  can  solve  equation  (2. 1)  effectively  when  ||®^/0y 


Is  small. 


A  generalization  of  equation  (3.1)  leads  fo  the  NLMS  methods  of  step  K 
In  the-  form 


VAh^‘VrhE  VAh>o„+r 

i=0 


(3.2) 


where  <*K  ^ °*  I  a0 1  + 1  ^Ko  (Ah))  |  >  °»  and  are  Independent  of  h,  but 

♦Ki  (Ah)  are  functions  of  h  and  nonsingular  A.  (The  Generalized  Adams- 

Bashforth  (GAB)  and  the  Generalized  Adams-Moulten  (GAM)  methods  for 
K  =  1,  2  are  suggested  by  the  work  of  Certalne.^) 

The  NLMS  methods  are  designed  to  solve  equation  (2)  effectively  when 
n8  V  -  i*-  >w\  is  large.  For  A,  which  is  a  constant  matrix,  we 
write  equation  (2. 2)  as 


£  (e"At  y)  -  e  At  g(t,y);  y(tQ)  =  yQ .  (3. 3) 


Integration  of  equation  (3.3)  over  the  Interval  ^tQ,  t^+,  J  gives 


m.  v  *Ah  ..  ,  ^  f  n+i  Afo+l-t')  j*, 

Y^W*®  y<tB)  +  i  ®  0(t’.y)dt'-  (3.4) 


n 
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Expressingfl/tr  . 


V  a°d,  next, 


Y(t 

n+-1 


iAb 


where 


We  define  the 


r'e  *V*£_^>do, 


n+i 

Xt<*h)=f  e  A(t  .,,, 

:  n+i  <"V*. 

n 

nonlinear  Multistep 


operator 

■n  fm  hj  „  r  „  Ah(K-i) 

^  <  yft  +  ih> 


v  K  -'i«  fywifcjtobe 


i  *  0 


(3.5) 


(3.6) 


^  *Ki  ^  9(t  +  ih,  y) 

Pw.  *  *  0 

Expanding 

ofh  att  =  t 


yields 


9(‘  +  ih.jr)  a  £ 

_  JI-*  (V  >*„»• 

(( we  substitute  »  a 


(3.7) 


j-0  -  ~  «••«> 

and  use  equation  (3.5)  for 
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-h 


Ah(K-i,elAhy  + 


t 


where 


C  (Ah),  f  .ieAh(K-'> 


i=0 


(Ah) 


J* 


K 


-h 


Z  wL  + 


j!  Ki 


(Ah) 


i=0 


A  consistent  NLMS  method  is  said  to  be  of  order  p  if 

CQ  =  Cj  =  •  •  •  =  Cp  =°.  but  Cp+1  • 

By  mathematical  induction,  it  is  seen  that 


x\  (Ah) 

7 


j2  =  o 


Therefore, 

-A<+1  Cj  (Ah,  e*h(K-‘»  £  + 

i=0  i«0 

•£  (l)J  *Kl  (Ah)  =0.  (3.9) 

i=0 
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A 


method  of  equation  (3.  2)  is  said  to  be  consistent 


3 


if 


K 


max 

n 


Ah  (K-i) 
e 


K 


Z>Ki<*h>«nJ 


is  small  as  h  —  0  . 


The  requirement  thatCj  (Ah)  =0  forj  K  0,  1,  .  .  p  yields  the  consistency 
and  permits  formulating  the  NLMS  method  in  the  following  matrix  form: 

14-  -  -HK4>  .  (3.10) 

where  E.  4- .  H  .  K,  and  $  are  described  by  the  expanded  forms  of  both 
explicit  and  implicit  forms. 


In  expanded  forms,  the  explicit  schemes  satisfy 


(3.11) 
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and  the  implicit  schemes  satisfy 


(3.12) 

We  determine  ♦  (Ah)  without  loss  of  generality  by  selecting  «  =  1 

Kl  K 


and  then  requiring  that  the  condition  of  strong  stability  be  realized  in  selecting 

a  ..  The  ♦  (Ah)  are  determined  using  the  above  matrix  formula,  which  can 
i  K.1 

be  considered  as  a  matrix  equation  for  the  K-step,  pth-order  method 

(K>1,  p  >0).  Below  we  show  that  ♦  (Ah)  can  be  determined  by.  first, 

Kl 

listing  the  formulas  for  the  remaining  +  (Ah),  and  second,  exhibiting  an 

Kl 

explicit  NLMS  method  of  order  2.  From  formula  (3.10)  for  nonsingular  H 
and  K  ,  we  get 


♦  =  -  K"1  H”1  !♦, 


(3.13) 


where  +  is  a  vector  of  dimension  (K-l)  or  K,  depending  upon  whether  the 
scheme  is  explicit  or  implicit.  Explicit  schemes  are  the  following: 

K  -  p  -  1.  ♦  (Ah)  =  -(Ah)"1  («0eAh  +  n  <3-14> 


9 
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K  =  p  =  2, 


♦  2  0(Ah,  =  -  (Ah,-2  [  ,#  (Ah-[)  e2**  -  a,  *h  -  a  ♦  Ahjp.,5, 
+2,i<Ah)  =  -(Ah)'2r«  e2Ah 

L  o 


+  a 


j  (I  +  Ah)eAh  +  (l+2Ah)]  (3. 


16) 


K  =  p  rr  3, 


*3 , =  "  <Ah)'3  [®  0  (I  -  |Ah+(Ah)2)e^h 


^  „  n  Ai,  2Ah 
+  <»1  fl  -  — )  e  +  «2 


Ah  .  Ah  .  3  *  2 1 

2  ^  6  +I  +  2^h  +  (Ah)  J  (3.17) 


♦  3  l<Ah)  =  -  (Ah)'3  [-2  *q  (I  -  Ah)e3Ah  +  a  1  (-21  +  (Ah)2)e2Ah 


-2  a 


2  "  +Al»  *h  -  PI  +  4Ah  +  3  (A I,)2)  ]  ,  p.,8, 

♦  3i2<*h)--<Ahf3  [«o(I-^-,e3Ah+«  p  +  *£*2 Ah 

1  Ci 

+  ff2  (I  +  fAh  +  Ah)2)eAh  +  (i  +|Ah  +  3  Ah)2)]  . 

Implicit  schemes  are  the  following: 

K«  P-  1, 

♦,  0  *h)  =  -  (Ahf2  [«o(Ah-I)eAh-i], 

*1.1  ,Ah)  =  '  <Ah)  2  [»„  e*h  *  p  +Ah)J  . 

K  =  p  =  2, 

♦*.0«*h>--l»h)'*  [•„  a  -f  Ah  ♦  (Ah,2,^ 

*2,1  (Ah)  “  '  [«Q  (-21  +  2Ah)e2Ah 

+  ai  <-2I+  <Ah)2)e  Ah  -  2  (I  +  Ah)]  , 


(3.19) 


(3 . 20) 

(3.21) 


(3.  22) 


(3.  23) 
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>2  2  (Ah)  =  -  <Ahf3  [  «0  fl  -  e24*1  +  „  (!+-&“,  e 

■1  z 

+  (1  +  jAh  +  (Ah)2)]  . 


Ah 


(3  •  24) 


K  =  p  =  3 , 

*3,0  <*h>  "  -  (*h)-4 


+al  (-I+Ah-J  <»' 

+  (-1  -  Ah  -  j  (Ah,2) 

♦  3  j  (Ah)  =  -  (Ah)’4 

+  (si-2Ah-i  (Ah)2  +  (Ah)  3j 
+  «2(3I+Ah-(Ah)2jeAh 

+  ^  31  +  4Ah  +  ~  (Ah)2 j 


'o  (-I  +  2Ah  -  ~  (Ah)2  +  (Ah)3  ) 


r2  (-1*! 

(Ah)2j  eAh 

(3 . 25) 

3Ah 

2Ah 

(3.26) 


*3,2  <*»'-<*!>) 


•„  (-31  +  4Ah  (Ah)2)  <..3Ah 

+  «J  (-31  +  Ah  +  (Ah)2)  e2**1 

+  *2  (•3I-2*'>  + j(Ah)2  +  (Ah)3)  eAh 
+  (-31  -  5  Ah  -  3  (Ah)2)l 


(3. 27) 


*3,3  <Ah)  '  -<Ah>'  [•0(I-Al,  +  5-«*h)2)e3Ah*a]  (I  -  J  (Ah)2)  e2**1 
+  *2(I+Ah+J  W»2KAh  +  (I  4  2  Ah  (Ah)2  +<Ah)3)]  (3. 28) 
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4.  STEP  SIZES  OF  NLMS  METHODS 

Nonlinear  multistep  methods  are  designed  to  avoid  the  use  of  small 

stepsizes where  g  (t,y)  Is  a  slowly  varying  function  that  can  be  approximated 

by  a  low-order  polynomial  in  t.  To  demonstrate  this,  we  give  an  analysis 

below  and  show  that  a  larger  step  size  can  be  selected  using  NLMS  methods 

than  using  LMS  methods.  From  equation  (3.  2),  since  a  ^  0,  we  can  write 

k  , 


n+k  a , 


K-l 

-  2>, 

i-0 


cAh<K~Vi+h  t  ***** 


which  is  of  the  form 


X=  G  (y), 


where  V  =  y  .  The  successive  iterative  form  gives 
n+k 

y(Y+1)  =  G  (y(V))  (4.2) 

for  any  initial  vector y^  . 

Let  G(y )  be  defined  for  ||  y  ||  <  *  ,  and  let  there  exist  a  constant  k 

such  that  0  <  k  <1.  Then,  G(y)  satisfies  the  condition 

II G  (  y*)  -  G  (y)||  <  k  ||  y  *-y||  .  (4.3) 

Using  the  definition  of  G  (y).  formula  (4.1),  and  the  fact  thatg(t,  y)  satisfies 

the  Lipschitz  condition  with  Lipschitz  constant  L,  we  see  that  condition  (4.3) 
is  satisfied  by 


h  II  *  KK 


(4.3) 


IS. 

for  sufficiently  small  h  and  for  all  II  A  II  <  «> 
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For  the  iterative  procedure  (4.2)  to  converge  for  arbitrary  initial 
y^,  k  is  required  to  be  less  than  1: 


k  <  1 


L  <  1. 


(4.5) 


Conventionally,  when  using  LMS  K-step  methods  with  ®K  =  1.  we 
select  h  such  that 


PK  hL*~  k  (<  1). 


(4.6) 


Similary,  for  NLMS  K-step  methods,  we  select  h^  to  satisfy  condition  (4.5): 


*KK  hN  L  k' 


(4.7) 


Combining  (4.6)  and  (4.7),  we  see  that 


♦kk  -1  (*hN) 


(4.8) 


For  ||  $KK  1  (Ah^  ||  not  too  small,  we  know  that  L*  »  L;  therefore, 
h^  >>  h.  This  tells  us  that  we  can  choose  a  much  larger  step  size  using  NLMS 
then  we  can  choose  using  LMS. 

The  quantity  h^  can  be  selected  to  satisfy 


advantage  is  demonstrated  by  the  problem  below. 


(4.9) 


This 


We  use  problem  1  from  section  8.3  to  demonstrate  the  step  size  choice. 
The  problem  is 

y'=  -iooy+  (1  +  t2);  y (0)  =  1  . 

Table  4-1  compares  two  methods  of  solution. 
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5.  FEATURES 


5. 1  VARIABLE  STEP  SIZE  AND  FIXED  STEP  SIZE 


Variable  step  size  Is  a  desirable  feature  that  Is  incorporated  into  this 
approach.  By  controlling  the  step  size,  the  solution  may  be  obtained  quicker 
through  the  use  of  a  larger  step  size.  In  some  instances,  the  user  may  desire 
to  examine  the  solution  at  a  specified  time  value.  The  variable-step-size 
technique  may  not  meet  this  requirement  because  a  larger  step  can  often  bypass 
that  point.  This  is  commonly  encountered  In  all  variable-step-size  programs. 
Although  NLMS  metl ods  are  designed  to  avoid  the  use  of  small  step  sizes, 
the  program  is  designed  with  this  option.  However,  the  user  can  still  use  the 
variable  step  size  if  he  desires  and  select  the  maximum  step  size  HMAX. 


To  exercise  this  option,  we  require  that  the  indicator  IPC  be  de¬ 
fined  as  follows: 


IPC  =  0:  fixed  step  size  with  INDEX  = 
IPC  /  0:  variable  step  size;  PCm. 


0  explicit 
1  implicit 


If  the  present  step  size  h  needs  a  change,  it  is  changed  by  an  amount 
rh  where  1/2  s  r  <  2.  To  save  computing  time,  we  change  h  by  halving 
or  doubling,  depending  upon  the  need.  This  is  handled  automatically  by  the 
program. 


An  example  to  demonstrate  this  feature  Is  given  in  problem  1  of  section 
8. 3.  Closely  connected  with  the  variable-step-size  procedure  is  the  PCm 
procedure,  which  will  be  described  in  the  next  section. 

5.  2  THE  PCm  PROCEDURE 


The  PC™  procedure  stands  for  "predict  and  correct  m  times."  This 
procedure  has  been  widely  used  to  achieve  efficiency.  Henrici2  discussed 
a  single  correction.  Others  feel  that  correction  may  be  needed  more  than 
once.  It  is  frequently  true  that  a  good  predictor  needs  but  a  single  correction. 
In  this  package,  we  set  the  upper  bound  of  m  to  be  3. 


It  should  be  noted  that  when  NLMS  methods  are  applied  to  solve  stiff 
equations,  when  9  is  a  function  of  t  alone  or  is  a  constant,  the  PCm  procedure 
may  not  be  needed. 
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5.3  SELF-STARTING 

The  possibility  of  missing  starting  values  is  common  with  all  multistep 
methods.  The  missing  values  have  to  be  provided  by  an  independent  method. 
The  Runge-Kutta  and  Adam-Bashforth  methods  are  often  used  as  starters. 

The  self-starting  procedure  is  provided  here  by  using  NLMS  or  LMS  first- 
order  methods.  As  mentioned  previously,  NLMS  methods  allow  the  use 
of  large  step  size  to  solve  stiff  equations;  therefore,  we  shall  adopt  the 
explicit  NLMS  methods  of  order  one  as  a  starter  for  solving  stiff  equations. 
The  Adams-Bashforth  first-order  method  is  used  as  a  starter  for  solving 
nonstiff  equations,  which  requires  the  use  of  a  very  small  step  size.  In 
the  above  case,  LMS  methods  are  called  for  and  the  self-starting  procedure 
requested.  It  is  suggested  that  the  user  specify  a  small  initial  step  size 
to  achieve  satisfactory  accuracy. 

In  most  instances,  the  use  of  NLMS  explicit  methods  of  order  one 
to  self-start  is  recommended  because  these  methods  are  a  good  starter  and, 
most  of  all,  because  the  step  size  is  not  restricted;  the  only  requirement 
is  that  the  g  function  must  be  defined.  To  do  this,  a  simple  modification 
can  be  Incorporated  in  the  START  subroutine.  The  user  can  save  the  indica¬ 
tor,  METHOD,  from  the  argument  and  set  it  to  1  after  entry,  then  reset 
METHOD  to  the  original  state  before  RETURN. 

5.4  SELECTION  OF  CHARACTERISTIC  COEFFICIENTS 

The  characteristic  polynomial  of  LMS  methods  takes  the  form 

K 

pu)=  E  Bi  ;i  =  0’aK  =  1-  (5*X) 

i=0 


Let 


be  the  loots  of  (5. 1). 


The  relationship 


K 


n 

i 


K 


i  =  o 


(5.2) 
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enables  us  to  express  as  a  function  of  4/.  the  4^  determine  the  root  condi¬ 
tion  and,  thus,  the  stability.  In  this  manner,  after  the  4^  are  chosen  to 
satisfy  the  root  condition  of  stability,  the  «  can  be  determined. 


For  NUN'  S  methods,  first  we  determine  (Ah),  which  are  dependent 

Ki 

on  a  as  given  by  formulas  (3. 14)  -  (3.  28).  After  $  (Ah)  are  determined, 
l  Ki 

the  numerical  solution  is  calculated  by  means  of  formula  (3.2),  which  is 

dependent  on  a  . 


For  LMS  methods,  the  a  are  usually  tabulated.  One  should  note 
that,  in  formula  (2.2),  if  A  =  O  ,  NLMS  methods  reduce  to  LMS  methods. 
If  we  multiply  (3 . 9)  by  [(Ah)J+1]"  and  let  ||  A  H  — 0,  we  get 


lJ  +  1 
(j  +1)1 


j! 


£ 


*’v»  • 


(5.3) 


For  j=0,  1,. . .  ,p-l,  equation  (5.3)  gives  a  formulation  of  the  LMS  methods 
of  order  p.  The  matrix  form  is 


Using  (5.4),  we  can  express  p .  (■♦  (Q)) 

*  Ki 

by  using  the  explicit  LMS  methods  of  order 


as  a  function  of  « 

2  whose  p  coefficients 


r 


tfj/2 

V2  +2*2 


)■ 


We  proceed 
are 


(5.5) 
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The  selection  of  two  roots,  0  and  1,  leads  to  strong  stability  and  gives 
a  0,  Ot  =  -1  and  a  -  1.  This  implies  that 

OX  ^ 


which  gives  the  Adams-Baahforth  method  of  order  two.  The  author  is  unaware 
of  any  program  that  allows  a  user  to  select  a  .  One  does  not 

1  1 

need  to  make  a  choice  when  using  tabulated  a  values.  He  has  found  that 

the  selection  of  roots  of  the  characteristic  polynomial  can  minimize  the 
initial  local  discretization  error  Cp+1  •  Therefore,  the  selection  of  a  }  may 

aid  in  studying  the  minimization  of  C  .  This  is  still  an  open  problem. 

p  +  1 

Hence,  this  option  may  not  be  immediately  exercised,  but  it  is  available  when 
needed . 


Different  characteristic  roots  satisfying  the  root  condition  give  the 
characteristic  coefficients  a  Different  families  of  implicit  NLMS  methods 

of  order  two  are  selected  to  solve  problem  3  in  section  8.  These  families, 
their  characteristic  roots,  and  the  associated  characteristic  coefficients  are 
identified  by  the  following: 

1 .  Strongly  Stable  Generalized-Adams  Family 
characteristic  roots:  0,1 

characteristic  coefficients:  a  Q  0,  a  ^  -  -1,  a  9  =  1. 

2.  Strongly  Stable  Family 
characteristic  roots:  1/2,  1 

characteristic  coefficients.-  a  1,  2,  O  -  -3/2,  a  ~  ~  1 

0  1  " 

3.  Weakly  Stable  Generali  zed -Milne  Family 
characteristic  roots:  -1,  1 

characteristic  coefficients:  C*  -1,  0=0,  0tn  =  1 

0  ’  1  2 
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5.5  A  AS  A  FUNCTION  OF  t 
Consider  the  initial  value  problem 

y'=  A(t)  y+  fl(t.y);  y(to)  =  yo  (5.6) 

over  the  interval  ft  ,  t  +  nh  1  .  We  periodically  decompose  equation  (5.6) 
into  L  0  0  J 

y'=  A(t.)  y  +  9(t.  y );  y(t 0)  =  yQ.  (5.7) 

where 

5(t,  y)  =  j  A(t)  -  A(t.)  j  y+g(t,y),  (5.8) 

such  that  Re  |  X  (A(t.))  |  <  0.  We  periodically  evaluate  A (t)  at  t  =  t..  The 
requirement  that  Re  j  X  (A(t.)) }  <  0  can  be  realised  by  construction  and  can 
be  assured  since  Re  I  X  (A(t))|  <  0  for  all  t.  The  property  of  g(t,  y)  can 
still  be  maintained  for  small  h  since  j|  A(t)  -  A(t})  ||  is  small  and  can  be  con¬ 
trolled  by  the  step  size.  In  cases  where  the  user  possesses  prior  knowledge 
about  the  problem,  he  can  introduce  a  constant  A  and  use  the  same  decomposi¬ 
tion  technique  and  perhaps  obtain  quicker  and  more  precise  results.  This  is 
not  to  say  decomposition  is  unique. 

5.6  INCLUSION  OF  LMS  METHODS 

As  a  part  of  this  package  LMS  methods  are  included,  whose  step 
numbers  are  chosen  to  be  compatible  with  the  same  step  numbers  used  for 
NLMS  methods  that  do  not  exceed  order  three.  The  application  of  LMS 
methods  is  indicated  by  an  indicator,  METHOD,  which  has  the  following 
definition: 


METHOD  = 


0,  LMS  methods 
1,  NLMS  methods. 


The  inclusion  of  LMS  methods  is  merely  an  option.  The  reader  is 
asked  to  consult  reference  4  for  an  efficient  implementation  of  the  Adams 
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methods.  However,  the  LMS  methods  incorporated  in  this  package  cover  a 
complete  linear  multistep  family  rathei  than  an  Adams  family  alone. 


5.7  MISCELLANEOUS  FEATURES 

1.  For  solving  stiff  equations,  ifg(t,y)  =  g  (t).  i.e.,  independent 
of  y,  no  correction  is  needed.  To  exercise  this  feature,  set  the  indicator 
EPC  equal  to  zero  so  that  unnecessary  computations  are  avoided. 


2.  If  the  stiff  equation  to  be  solved  is  homogeneous  and  A  is  a  con- 

K 

stant  matrix,  there  is  no  need  to  calculate  h  7  $,..(Ah)f|  By  setting  the 

Ki  »n+i 


i=0 


K 


indicator  IGFN  equal  to  zero,  the  program  will  bypass  the  h  X^Ki^fln+i 

1=0 


computation  loop.  In  the  event  that  A  is  a  function  of  t,  even  if  the 

system  is  homogeneous,  theng(t,  y)  =  0  must  be  defined  in  the  GFN  subroutine. 

3.  A  user's  starter  can  be  supplied  to  replace  the  START  through 
either  the  argument  or  the  complete  replacement  of  the  subroutine. 

Subroutine  INVERT  can  be  treated  in  the  same  manner  as  above. 
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6.  PROGRAMMING  ASPECTS 


6. 1  NOTES  TO  USERS 


1.  Double  precision  arithmetic  is  used  for  all  computations. 

2.  Subroutines  not  used  by  the  user  can  be  Ignored.  They  are  well- 
defined  Internally. 

3.  Subroutines  that  are  required  to  be  supplied  oy  the  user  must  no* 
be  ignored  in  order  to  obtain  correct  computations. 

a.  If  NLMS  methods  are  used,  GFN  must  be  supplied,  except 
fl(t,y)  =0(IGFN  =  0). 

b.  If  LMS  methods  are  used,  FFN  must  be  supplied. 

c.  If  A  is  a  function  of  t,  AFNT  must  be  defined. 

4.  The  solution  vector  is  either  in  YNEW(I)  or  Y(KSTEP  +  1,  i).  The 
associated  tj  is  in  TEA.  The  user  who  desires  to  output  his  results  can 
incorporate  the  desired  output  format  into  PRINT  subroutine. 

5.  The  existing  Univac  1108  EXEC  8  library  matrix  inversion  sub¬ 
routine  DGJR  is  used  for  test  problems. 

6.  Certain  detections  to  inputs  are  provided  in  FORTRAN  V  programs. 
Default  values  are  assigned  to  protect  a  great  number  of  Inputs.  This  advan¬ 
tage  can  be  used  since  the  inputs  are  defined  on  NAMELIST.  However,  this 
advantage  is  not  applicable  to  ANSI  FORTRAN  users;  their  inputs  must  be 
supplied  correctly. 


7.  Meanings  of  a  few  special  indicators  in  DIFEQ: 


IG 

IV 

INDEX 

IAT 


Indicators,  brought  to  the  DIFEQ  program  through  the 
CALL  argument,  are  IG,  IV,  INDEX,  IAT.  These 
are  normal  Indicators.  I 


indicates  whether  or 
not 


g(t,  y)  =0  . 

a  variable  step  size  is  used, 
a  method  Is  explicit. 

A  is  a  function  of  t. 
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Other  than  these  indicators,  there  are  a  few  indicators  specially  used  by 
DIFEQ  to  produce  computational  efficiency  and  accuracy.  Their  meanings 
need  to  be  explained. 

LMT.  Normal  use  is  to  define  LMT  =  m  ~  3  when  a  variable-step-size, 
m  - 

PC  procedure  is  used.  To  ensure  this,  LMT  is  set  to  3  initially.  In  the 
event  a  fixed-step-size,  implicit  method  is  asked  for,  we  make  use  of  this 
indicator  to  yield  the  following  control. 

a.  When  IV  0  and  INDEX  =  1,  a  fixed-step-size,  implicit  method 
is  asked  for. 

b.  After  a  complete  calculation  of  the  required  implicit  formula, 
LMT  is  changed  to  contain  1.  The  purpose  of  this  modifica¬ 
tion  is  to  inform  the  program  to  predict  the  next  yn  value 

in  order  to  implicitly  calculate  theyn  .  This  internal  modifica¬ 
tion  does  not  affect  the  subsequent  calculations. 

IT  EH.  A  counter  for  every  successful  corrector's  convergence  or 
every  fixed-step  calculation.  To  take  advantage  of  this,  this  indicator  is 
designed  to  appear  in  the  argument  of  PRINT,  so  that  the  user  may  use  this 
information  to  control  his  printouts.  As  an  example,  if  the  user  wants  to 
print  out  results  at  every  100  successful  calculations,  he  npeds  to  define  a 
statement  in  the  PRINT  program  that  reads: 

IF(MOI)(ITER,  100)  .  EQ.  0)WRITE(4, 12)H,T,  (Y(L),  L=1,N). 

ISTEP.  If  a  fixed  step  size  is  used  to  proceed  for  subsequent  computa¬ 
tions,  the  matrix  exponentials  and  the  determination  of  $  (Ah)can  be 

Ki 

calculated  only  once.  When  a  variable  step  size  is  used,  it  is  necessary  to 
calculate  the  matrix  exponentials  and  to  determine  the  4  (Ah)  for  every  h. 

When  the  program  finds  a  successful  h  that  leads  to  the  corrector's  convergence 
as  well  as  satisfying  user’s  tolerance,  we  take  a  new  approach:  We  can  use  the 
same  h  to  proceed,  if  required  conditions  are  met,  for  at  most  three  calcula¬ 
tions;  we  then  consider  doubling  the  step  size.  Suppose  we  double  the  step 
size  too  soon.  Then  it  is  necessary  to  calculate  the  matrix  exponentials  and 
to  determine  the ♦  (Ah).  If  this  new  h  does  not  lead  to  the  corrector's 

convergence,  all  above  computations  are  wasteful.  Then  we  need  somehow 

to  use  the  matrix  exponentials  and  the  4  (Ah) for  the  old  h.  Without 

Ki 

worrying  about  memory  space,  the  old  values  can  be  saved;  without  worrying 
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about  computation  time,  the  old  values  can  be  recalculated.  Certainly,  we 
cannot  save  both  the  time  and  the  space  simultaneously.  Now,  the  reason 
that  ISTEP  is  used  to  minimize  cost  becomes  apparent.  In  the  case  of  solv¬ 
ing  stiff  equations,  chances  are  good  for  success,  even  if  h  is  doubled. 

The  user  can  take  advantage  of  this  by  requiring  the  limit  of  ISTEP  to  be  1. 

The  change  to  the  program  Is  almost  trivial.  The  present  setup  doubles  h  after 
each  successful  calculation. 

IMIN.  To  perform  floating-point  addition  between  T(l)  and  H,  the 
difference  between  T(l)  and  H  must  be  within  the  allowable  significant 
digits  of  the  computer.  To  maintain  m  significant  digits  in  floating-point 
addition,  the  exponent  (T(l))  minus  the  exponent  (H)  cannot  exceed  m.  We 
define  HMIN  to  mean  h  and  to  satisfy  this  range.  In  reality,  even  if  the 

Lipschitz  constant  is  large  and  even  if  the  LMS  methods  are  called  for, 

HMIN  will  not  fall  outside  of  the  significant  digit  range,  and  we  can  assign 
a  larger  h  ^  to  avoid  excessive  computations.  Once  this  safeguard  is 
built  in,  wffen  H  reaches  HMIN  more  than  three  times,  we  will  not  continue. 

We  will  terminate  the  program  with  a  message  for  the  user  to  tell  him  what  to 
do.  The  variable  IMIN  Is  defined  to  count  how  many  times  H  reached  HMIN. 

8.  Even  though  programs  work  for  the  same  computer  in  different 
locations,  the  Input/output  (I/O)  units  may  differ  In  number;  users  are  ad¬ 
vised  to  make  sure  that  the  unit  numbers  are  in  agreement  with  the  program. 

6.  2  FORMATS  OF  USER-SUPPLIED  SUBROUTINES 

A  few  optional  subroutines  must  be  supplied  by  the  user  under  certain 
conditions.  These  conditions  are  listed  below  along  with  their  required 
formats. 

1.  INVERT.  This  is  a  matrix  inversion  subroutine,  called  by  NLMS 
methods  whose  CALL  statement  has  the  format 

CALL  INVERT  (A,  N,  B), 


where 

A  is  an  (n  x  n)  matrix 
N  is  the  order  of  matrix  A 
B  contains  the  A  inverse  after  exit. 


25 


TR  5341 


As  an  example,  using  an  existing  Univac  1108  library  subroutine  DGJR  (Gauss- 
Jordan  Reduction),  the  structure  of  INVERT  should  look  like 


SUBROUTINE  INVERT  (A,  N,  ANS) 
PARAMETER  NR=20,  NC=20 
IMPLICIT  REAL*8  (A-H,  O-Z) 

DIMENSION  A  (NR,  NC),  ANS  (NR,  NC) 
DIMENSION  V  (2),  JC  (NC) 

V  (1)  =  1  .DO 
DO  1  I  =  1,  N 
DO  1  J  =  1,  N 

1  ANS  (I,  J)  =  A  (I,  J) 

CALL  DGJR  (ANS,  NR,  NC,  N,  N,  $10,  JC,  V) 
RETURN 
10  WRITE  (4,  2) 

2  FORMAT  (3X,  'MATRIX  INVERSION  ERROR’) 
RETURN 

END 


2.  BEGIN.  This  subroutine  supplies  initial  values.  It  is  required  if 
the  self -starting  procedure  is  not  used.  Its  CALL  statement  has  the  format 

FORTRAN  V:  CALL  BEGIN  (K,  H,  N,  Y,  YN,  YOLD,  INVERT, 

IT,  METHOD,  IAT) 

ANSI  FORTRAN:  CALL  START  (K,  H,  N,  Y,  YN,  YOLD,  IT, 
METHOD,  IAT), 


where 


K  is  the  step  number 
H  is  the  step  size 
N  is  the  order  of  the  system 

Y  is  a  two-dimensional  array  containing  the  initial  values;  at  most, 
4  rows,  20  columns. 

YOLD  is  the  same  dimension  as  Y,  used  as  a  working  storage. 
INVERT  is  the  name  of  a  matrix  inversion  subroutine. 
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IT  is  an  indicator  that  indicates  whether  the  calculation  of  a 
g  -function  is  needed. 


METHOD  indicates  use  of  either  NLMS  or  LMS  explicit  of  order  one. 


IAT  indicates  whether  A  is  a  function  of  t. 


As  an  example,  if  one  uses  exact  initial  values  for  problem  5,  section  8,  the 
structure  of  BEGIN  should  look  like  this  in  FORTRAN  V: 

SUBROUTINE  BEGIN  (K,H,N,  Y,  YN.YOLD, INVERT,  IT,  METHOD, IAT) 
PARAMETER  NR  =  20,  NC  =  ?0 
IMPLICIT  REAL*8  (A=H,  O-Z) 

EXTERNAL  INVERT 

COMMON  A  (NR,  NC),  ALPHA  (4),  T  (NC) 

DIMENSION  Y  (4,  NC),  YN  (N),  YOLD  (4,  NC) 

H2  *  T  (1)  +  H 
H3  =  H2  +  H 

Y(2, 1)  =  DEXP  (-H2*H2/2.)  -  DEXP  (-H2)  +1.0 
Y(3 , 1)  =  DEXP  (-H3+H3/2.)  -  DEXP  (-H3)  +  1.0 
RETURN 
END. 

In  ANSI  FORTRAN,  the  structure  should  be  as  follows: 

SUBROUTINE  START  (K,H,N,Y,  YN.YOLD, IT,  METHOD, IAT) 

COMMON  A  (20,  20),  ALPHA  (4),  T  (20) 

DIMENSION  Y  (4,20),  YN  (20),  YOLD  (4,20) 

DOUBLE  PRECISION  A,  ALPHA,  H,  T,  Y,  YN,  YOLD,  H2,  H3 
H2  =  T  (1)  +H 
H3  =  H2  +  H 

Y  (2,  1)  =  DEXP  (-H2*H2/2.D0)  -DEXP  (-H2)  +1.D0 

Y  (3,1)  -  DEXP  (-H3*H3/2.D0)  -DEXP  (-H3)  +1.D0 
RETURN 

END 


3.  GFN.  This  subroutine  calculates  the  g(t,  y )  at  t  «  t  for  every 

i.  It  is  required  by  NLMS  methods  only  if  the  differential  equation  is  non- 
homogeneous.  Its  CALL  statement  has  the  following  format: 

FORTRAN  V:  CALL  GFN  (G,  H,  N,  Y,  K,  T,  A,  NR,  NC) 

ANSI  FORTRAN:  CALL  GFN  (G,  H,  N,  Y,  K,  T,  A), 
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where 


G  is  a  vector  that  contains  calculated  g  function  values 
H  equals  the  step  size 
N  is  the  order  of  the  system 

Y  is  a  two-dimensional  array  Y  (I,  J):l  ^  I  <  4; :  l^J^N 

K  equals  the  step  number 

T  is  a  vector  that  contains  time  values 

A  is  an  (n  x  n)  matrix 

NR  equals  the  number  of  rows  of  A 

NC  is  the  number  of  columns  of  A. 

For  example,  to  solve  a  system  of  two  stiff  equations,  problem  3  of  section  8, 
the  structure  of  the  GFN  should  be  as  follows  in  FORTRAN  V: 

SUBROUTINE  GFN  (G,  H,  N,  Y,K,  T,  A,  NR,NC) 

IMPLICIT  REAL*8  (A-H,  O-Z) 

DIMENSION  A  (NR,  NC),  Y  (4,  NC),  G  (NC),  T  (NC) 

G  (1)  =  1.0 
G  (2)  =  1.0 
RETURN 
END 

The  structure  in  ANSI  FORTRAN  should  be  as  follows: 

SUBROUTINE  GFN  (G,H,N, Y,K,T,  A) 

DOUBLE  PRECISION  A  (20,  20),  Y  (4,  20),  G  (20),  T  (20),  II 

G  (1)  a  1  .DO 

G  (2)  =  1.D0 

RETURN 

END 
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4.  FFN.  This  FFN  subroutine  calculates  f  (t,  y),  which  is  required 
only  byLMS  methods.  Its  CALL  statement  has  the  format 

CALL  FFN  (Y,N, FN, I), 


where 


Y  is  a  2-dimensional  array  Y  (I,  J); 

1  *  I  <  step  number  +  1  ;  1  <  J  <  N 
N  equals  the  order  of  the  system 
FN  is  a  vector  to  store  calculated  f  function  values 


I  is  the  intermediate  step  index  . 


For  example,  to  solve  a  system  of  two  equations,  problem  4  of  section  8, 
using  LMS  methods,  the  structure  of  FFN  should  be  as  follows  for 
FORTRAN  V: 

SUBROUTINE  FFN  (Y.N.FN.I) 

PARAMETER  NR  =  20,  NC  =20 
IMPLICIT  REAL*8  (A-H,  O-Z) 

COMMON  A  (NR,  NC),  ALPHA  (4),  T  (NC) 

DIMENSION  FN  (1),  Y  (4,  NC) 

FN  (1)  =  Y  (I,  2) 

FN  (2)  =  -Y  (I,  1) 

RETURN 

END 

For  ANSI  FORTRAN,  the  structure  is 

SUBROUTINE  FFN  (Y.N.FN.I) 

COMMON  A  (20,  20),  ALPHA  (4),  T  (20) 

DIMENSION  FN  (20),  Y  (4,  20) 

DOUBLE  PRECISION  A,  ALPHA,  FN,  T,  Y 
FN  (1)  =  Y  (I,  2) 

FN  (2)  =  -Y  (I,  1) 

RETURN 

END 
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5.  AFNT.  This  AFNT  subroutine  should  be  supplied  by  the  user 
when  A  is  a  function  of  t.  This  subroutine  calculates  every  A  (t)  at  t  -  t  . 

Its  CALL  statement  has  the  format 

FORTRAN  V:  CALL  AFNT  (A,N,T,NR,NC) 

ANSI  FORTRAN:  CALL  AFNT  (A,N,T), 


where 


A  is  a  two-dimensional  array,  A  (NR,  NC) 

N  equals  the  order  of  the  matrix 
T  is  a  time  variable 

NR  is  the  number  of  maximum  rows  of  A 

NC  equals  the  number  of  maximum  columns  of  A  . 

For  example,  to  solve  a  system  of  one  equation  (problem  5)  using  an  NLMS 
method  of  order  2,  the  structure  of  AFNT  in  FORTRAN  V  should  be: 

SUBROUTINE  AFNT  (A,  N,  T,  NR.NC) 

DOUBLE  PRECISION  A  (NR,  NC),  T 

A  (1,  1)  =  “T 

RETURN 

END 

In  ANSI  FORTRAN,  it  should  be: 

SUBROUTINE  AFNT  (A,N,T) 

DOUBLE  PRECISION  A  (20,  20),  T 

A  (1,1)-  -T 

RETURN 

END 


6.  PRINT.  This  PRINT  is  designed  for  the  user’s  convenience  to 
output  the  results  in  his  desired  formats.  Its  CALL  statement  has  the  format 

CALL  PRINT  (H,T,Y,N,I), 


i 

i 

i 

t 

1 


1 
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where 


H  is  the  present  step  size 
T  is  the  present  time  step 

Y  is  a  one-dimensional  array  containing  the  current  solution 
N  is  the  order  of  the  system 
I  is  an  iteration  counter. 

For  example,  results  are  required  to  be  printed  out  at  t  =  5  and  10  of  problem 
5  (section  8)  along  with  exact  solutions  at  the  same  time  values.  The  PRINT 
statement  can  be  designed  as  follows: 

SUBROUTINE  PRINT  (H,T,Y,N,I) 

DOUBLE  PRECISION  H,  T,  Y  (1) 

12  FORMAT  (IX,  2  (E15.8,  2X),  4X,  5  (E15.8,  2X)) 

IF  (T-4,  995D0)  84,84,83 

83  IF  (T-5.  005D0)  1184,1184,85 

85  IF  (T-9.  995 DO)  84,84,86 

86  IF  (T-10.  005D0)  1184,1184,84 
1184  WRITE  (4,12)  H,T,  (Y(J),  J  =  1,  N) 

CALL  EXACT  (T.Y) 

84  RETURN 
END 
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7.  USER'S  GUIDE 

The  following  table  gives  the  options,  the  default  conditions,  and  the 
input  indicators  available. 

Table  7-1.  User's  Guide 


Option 

Default 

Condition 

Input  Indicator 

A  is  a  function  of  t 

0 

IAT  4  0 

ERR 

. 

0.1D-11 

F -function  (FFN) 

0 

Called  by  LMS  only  if 
METHOD  =  0 

final  time 

TMAX 

I  PC  = 

1  and 

fixed  step  size 

INDEX  = 

0  Explicit 

1  Implict 

G-function  (GFN) 

0 

=  0  not  needed 

IGFN 

%  0  need  for  stiff 
eqs. 

INDEX 

0 

as 

0  Explict 

1  Implitict 

initial  time 

T  (1) 

INVERT 

user-supplied  matrix 
inversion 
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Table  7-1.  User's  Guide  (cont) 


Option 

Default 

Condition 

Input  Indicator 

KSTEP,  step  number 

1 

KSTEP 

METHOD 

1 

0  LMS 

1  NLMS 

Pcm 

0 

IPC  =  1  or  IAT  4  0, 
m  =  3 

START 

self-starting 

User-supplied  starter 

step  size 

0.01 

H 

step  size  (maximum) 

0.1 

HMAX 

variable  step  size 

IPC  =  1 

initial  vector 

YZERO 

«i 

GAB-GAM 

ALPHA 

order  of  the  system 

N 
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8.  NUMERICAL  RESULTS  AND  INPUT  SETUPS 


8. 1  PROBLEM  IDENTIFICATION 

An  ordinary  differential  equation  is  said  to  belong  to  the  class  (N,  S, 
H,  L)  if  it  satisfies  the  following  definitions: 


N  * 


S  = 


0  first  order 
1  higher  order 
0  nonstiff 
1  stiff 


H  = 


0  homogeneous 
1  nonhomogeneous 


L  * 


0  linear  in 
1  nonlinear  in . 


For  example,  the  equation 


is  a  first-order,  nonstiff,  homogeneous,  linear  differential  equation.  There 
fore  it  belongs  to  (0, 0, 0, 0). 


8.2  SELECTION  OF  TEST  PROBLEMS 

6-11 

Seven  test  problems  have  been  selected  from  various  sources 
in  order  to  illustrate  the  different  options.  At  least  one  feature  is  exerc  se 
in  each  problem.  The  exactstart  will  be  used  for  most  test  P^btomB.  T 
formats  of  user-supplied  subroutines  are  described  within  each  problem 
FORTRAN  V  langua^.  The  main  body  of  ANSI 
FORTRAN  V  of  each  subroutine;  therefore,  the  ANSI  FORTR 
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omitted.  For  exact  ANSI  FORTRAN  formats,  the  user  can  consult  the  ANSI 
FORTRAN  in  the  appendix.  Table  8-1  itemizes  the  seven  problems, 
which  are  then  each  discussed  in  detail. 


8.3  PROBLEMS,  INPUT  SETUPS,  AND  SOLUTIONS 
Problem  1 


y*=  -loo  y+  (i  +  t  );  y(o)  =  1. 


Exact  Solution 


y<t»  =  <1-^-^>e',00t  +  ^+T^(10°2 12 

-200t  +  2). 


Classification 

(0,1, 1,0). 

Eigenvalues  of  A 

(-100  } 

Method  Applied 

Implicit  NLMS  of  order  3  with  fixed  step  size.  Linear  multistep  of 
order  3  with  PCm  procedure. 

Step  Size  Used 

1.25  and  2.5  for  NLMS.  Variable  slep  size  for  LMS  with  initial 
h  =  0.01. 

Time  Interval 

[0,  10] 

Special  Features  Exercised 

Fixed  step  size. 

Variable  step  size. 
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Special  Features  Exercised  (cont) 

Exact  start. 

PCm  procedure. 

FORTRAN  V  Inputs  Set-Up 

Fixed  step  size,  implicit  NLMS  of  order  3. 

$INPUTS 

N=  1,  KSTEP=3,  ALPHA(1)*O.DO,  O.DO,  -1.D0,  1.D0, 
T(1)*0.D0,  TMAX=.  1D2, 

YZERO(l)=l  .DO,  A(1,1)=-1.D2, 

H=1.25DO, 

IGFN=1,  IPOO,  METHOD=  1,  INDEX-1, 

$END 

Variable  step  size,  LMS  of  order  3. 

$ENPUTS 

N=  1 ,  KSTEP=3,  ALPHA(1)=0.D0,  O.DO,  -1.D0,  1.D0, 
T(l)=  ODO,  TMAX= .  1D2 , 

YZERO(1)=1.DO,  ERR= .  ID  -11, 

Hs.lD-3,  HMAX= .  ID-2 , 

IPC=1,  METHODS,  INDEX=0 
$END 

User- Supplied  Subroutines 


SUBROUTINE  GFN  (G,H,N,Y, J, T,  A,NR, NC) 
IMPLICrr  REAL*8  (A-H,  O-Z) 

DIMENSION  A  (NR,  NC),  Y  (4,  NC),  G  (NC),  T  (NC) 

TH  e  T  (1)  +  (J-l)  *H 

G  (1)  -  1 .  +  TH  *  TH 

RETURN 

END 
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User-Supplied  Subroutines  (co nt) 


SUBROUTINE  FFN  (Y,N,FN,I) 
PARAMETER  NR=20,  NC=20 
IMPLICIT  REAL*8  (A-H,  O-Z) 

COMMON  A  (NR,  NC),  ALPHA  (4),  T  (NC) 
DIMENSION  FN  (1),  Y  (4,  NC) 

FN  (1)  «  -100,  *Y  (I,  1)  +  (l.+T  (I)  **2) 

RETURN 

END 


An  Exact  Starter  -  BEGIN 


SUBROUTINE  BEGIN  (K,H,N,  Y,  YN.YOLD,  INVERT,  IT,  METHOD,  I  AT) 
PARAMETER  NR=20,  NC=20 
IMPLICIT  REAL*8  (A-H,  O-Z) 

EXTERNAL  INVERT 

COMMON  Y  (4,  NC),  YN  (N),  YOLD  (4,  NC) 

H2  =  T  (1)  +  H 

H3  *  H2  +  H 

H4  *  H3  +  H 

X  =  1.  -.  01-2.  /100.**3 

X4=X*DEXP(-100.*H4) 

X3*=X*DEXP(-100.*H3) 

X=X*DEXP(-100.  *H2) 

W=(100.  *H2)**2-200.  *H2  +  2. 

Y(2 , 1)=X  +  .  01+W/100.**3 
W=  (100.  *H3)**2-200.  *H3+2. 

Y(3 , 1)=X3+  .  01 +W/100.  **3 
W«(100.*H4)**2-200.*H4+2. 

Y(4 , 1)«X4+ .  01+W/100.  **3 

RETURN 

END 


SOURCE:  Lee1. 
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Numerical  Results 


Give  results  around  t  =  5,  10. 
Implicit  NLMS  methods  of  order  3. 


t 


y(t) 


.50000000+01 


.25900200-00* 
. 25900200 -00f 


.10000000+02 


LMS  methods  of  order  3 


. 10080020+01* 
.10080020+01 t 


t 

.50047500+01 

,10020750+02 


y(t) 

.25497628-00* 

. 25497628-OOf 

.10121522+01* 

.10121522+01 t 


Remarks 


Nonlinear  multistep  methods  can  arrive  at  t  **  10  in  one  step  by  using 
h  =  2. 5  at  the  cost  of  4. 1  seconds  CPU  time.  However,  in  order  to  obtain 
y  (5),  we  use  h  *  1.25;  the  CPU  time  costs  a  little  more— 4.4  seconds. 

These  NLMS  methods  are  designed  to  allow  the  use  of  large  step  size  and  to  be 
effective  for  those  g(t,  y )  belonging  to  the  class  of  low-order  polynomials. 
Therefore,  it  is  not  surprising  that  the  implicit  NLMS  method  of  order  3 
produces  accurate  results. 


♦Solution  produced  by  the  methods 
t  Exact  solution 
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Problem  2 


using  6  k  =  1. 0075 

p  =  p  =  0.0075 


v  =  0.075 

K\ 


J|  =  10 


-6 


Exact  Solution  T 

y  (2)  =  (-.6414  1073-07,  -.8552  1424  +  00) 

y(5)  *  (-.5130  4190-07,  -.6840  5581  +  00 ) 
y(10)_  ( -.3536  0130-07,  -.4714  6837  +00  ) 

r.lassificiation 

(0, 1,0,0) 

Eigenvalues  of  A 

{-.0744375,  -106  } 


Method  Applied. 

Implicit  NLMS  of  order  2  with  fixed  step  size. 

step  Size  Used  _ 

1.0 

Time  Interval 
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Special  Features  Exercised 

Fixed  step  size 
Exact  Solution 

FORTRAN  V  Inputs  Setup 

Fixed  step  size,  implicit  NLMS  of  order  2 


$INPUTS 

N=  2,  KSTEP  =  2,  ALPHA  (1)  =  O.DO,  -1.D0,  1.D0 
T  (1)  =  O.Du,  TMAX  -  10.  DO 
YZERO  (1)  =  1  .DO,  -1. DO,  H=1.D0 
IGFN  =  0,  IPC  =  0,  METHOD  -  1,  INDEX  =  1 
$LND 

User-Supplied  Subroutine 

Since  g(t,  y  )  =  O  .  IGFN  is  set  to  0.  Hence,  the  GFN  is  not  called  for. 
Source  Weaver10 


The  infinite-medium  reactor  kinetic  equations,  in  standard  form ,  with 
m  delayed  neutron  groups,  can  be  expressed  as 

m 


dn 
dt  " 

6  k  -P 
i 

■ n+  £  v. 

i  =  0 

dC 

i 

Pi 

"  ■  \  cf 

dt 

“  i 

where 

{Jj  =  the  fraction  of  the  total  number  of  fission  neutrons  belonging  to  the  i^ 
group  precursor  that  are  delayed , 

P  =  the  fraction  of  the  fission  neutrons  that  are  delayed, 

<L  =  neutron  generation  time  (set  ). 
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X. .  =  decay  constant  for  the  i  group  precursor, 
n  =  neutron  density, 

C.  =  the  concentration  of  the  i^  precursor,  and 
6^  =  a  constant  for  a  reactor  with  a  constant  excess  of  reactivity. 
Considering  a  single  delayed  neutron  group,  we  obtain 


Numerical  Results 


Give  results  at  t  =  2,  5,  10. 


t 

Y1  (t) 

y2(t) 

.2000  0000  +  01 

-.6414  1073-07 
-.6414  1073-07 

-.8552  1424+00* 
-.8552  1424+00 t 

.5000  0000+01 

-.5130  4190-07 

-.5130  4190-07 

-.6840  5581+00 
-.6840  5581+00 

.1000  0000+02 

-.353  6  0130-07 
-.3536  0130-07 

-.4714  6837  +  00 
-.4714  6837  +  00 

Remarks 


A-stable  NLMS  methods  can  solve  homogeneous  stiff  equations  exactly 
in  the  absence  of  round-off  errors.  An  extremely  small  step  size  is  required 
for  the  same  accuracy  {.lx  10-6)  if  LMS  methods  are  used. 


♦Produced  by  Implicit  NLMS  methods  of  order  2 
t  Exact  solution 
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Problem  3 


0  1 


V  =  V  10  -9 


F.vaet  Solution 


y(t)  = 


2e*-l 
2e^_l 


Classification 

(0,  0,  1,  0) 

Eigenvalues  of  A 

|-io,  l| 

Methods  Applied, 

Explicit  NLMS  method  of  order  3. 

step  Size  Used 

0.1  and  1.0. 

Time  Interval 

[0,  10] 


jcial  Featu™H  Exercised. 
Fixed  step  size. 
Self-start. 
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FORTRAN  V  Inputs  Set-Up 

Fixed  step  size,  explicit  NLMS  of  order  3. 


$INPUTS 

N=2,  KSTF,P«=3,  ALPHA  <1)=0. DO,  0.D0,  -1.D0,  1.D0, 
T(1)=0 . DO,  TMAX= .  1D2 , 

YZERO(l)=  1 .  DO,  1.D0, 

A(1 ,  1)  =  0 .  DO ,  A(1,2)=1.D0,  A(2, 1)=  .  1D2 ,  A(2, 2)= -9.DO, 
H*l.DO, 

IGFN=  1 ,  IPC=0,  METHOD- 1 ,  INDEX=0 
$END 

User-Supplied  Subroutines 


SUBROUTINE  GFN  (G,H,N,Y,  J.T.A.NR.NC) 
IMPLICIT  REAL+8  (A-H,  O-Z) 

DIMENSION  A  (NR,  NC) ,  Y  (4,  NC),  G  (NC),  T  (NC) 

G  (1)  =  1.0 

G  (2)  =  1.0 

RETURN 

END 


Source:  Lee  . 


Numerical  Results 


Gives  results  at  t=l,  5,  and  10. 


t 

.10000000+01 

.50000000+01 

.10000000+02 


Yj(t) 

.44365637+01 

.44365637+01 

.29582632+03 

.29582632+03 

.44051932+05 
.44051932  +05 


y2(t) 

.44365637+01* 
.443  65637+ 01  f 

.29582632+03 

.29582632+03 

.44051932+05 

.44051932+05 


♦Produced  by  explict  NLMS  methods  of  order  3. 
t Exact  solution. 
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Since  LMS  methods  are  not  used,  FFN  is  not  called  for.  The  built-in 
self-starter  is  used. 


Problem  4 

u"  +  u  =  0;  u  (0)  =  0,  u *  (0)  =  1. 
Equivalent  Problem 


Exact  Solution 


u  =  sin  t. 


Classification 


Eigenvalues  of  A 


(1, 0,0,0) 

l*i  I  • 


Method  Applied 

LMS  methods  of  order  2  with  PCm  procedure. 
Step  Size  Used 

Initial  step  size  ■  *  . 


Remarks 

In  the  absence  of  round-off  errors,  NLMS  methods  of  different  orders 
are  equivalently  accurate  if  g(t,  y>  is  a  constant;  NLMS  methods  of  any  order 
should  produce  accurate  results  for  this  problem. 
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[O,  6  it]  . 


Special  Features  Exercised 
Variable  step  size 
PC  procedure 


Exact  start. 

FORTRAN  V  Inputs  Set-Up 
$INPUTS 

N=2,  KSTEP=2 ,  IPC*1,  INDEX* 0,  METHOD* 0 
T(1)*0.D0,  TMAX*  18.7 5D0 
YZERO(1)=0 .DO,  1.D0,  ERR=.1D-11 
ALPHA(l)  =  0. DO,  -1. DO,  1.D0 
$END 

In  MAIN,  H  is  defined  to  be  tt  initially  and  HMAX  *  H/16. 


User-Supplied  Subroutines 


SUBROUTINE  FFN  (Y,N, FN,  I) 
COMMON  A  (20,  20),  ALPHA  (4),  T  (20) 
DIMENSION  FN(20,  Y(4,20) 

FN  (1)  =  Y  (I,  2) 

FN  (2)  »  -Y  (I,  1) 

RETURN 

END 


Source: 


Krough 


9 


Numerical  Results 


Give  maximum  absolute  error  at  t  ■  2,  4,  6. 

t  Maximum  Absolute  Error 

2  w  .2307  1215-12 

4  w  .2307  0955-12 

6  w  .2307  1128-12 
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Remarks 


This  problem  demonstrates  that  high-order  differential  equations  can 
be  decomposed  into  a  system  of  first-order  equations  so  that  this  package  is 
applicable.  Since  this  is  not  a  stiff  equation,  we  choose  the  variable -step-size 
and  PCm  procedure  to  examine  the  efficiency  of  the  LMS  portion. 


Problem  5 

y'=  t  (1  -  y)  +  (1  — t)e_t;  y  (.1)  =  1.0901751. 


Equivalent  Problem 

y'=  -t  y  +  t  +  (l-t)e-t;  y(.l)  =  1.0901751. 


Exact  Solution 


Classification 


0  <  t  <  1  (0,0, 1,0) 

1<  t  (0,1,1, 0). 


Eigenvalues  of  A 


I 


Methods  Applied 

NLMS  methods  of  order  3. 

Step  Size  Used 

Variable  step  size  with  initial  h  =  0.01. 
Time  Interval 

[0.1,  50. 0]  . 
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Special  Features  Exercised 


Automatic  periodic  decomposition  of  A(t). 
PCm  procedure. 

Variable-step-size  procedure. 

FORTRAN  V  Inputs  Set-Up 
Variable  step  size,  NLMS  of  order  3 


$  INPUTS 

N=s  1 ,  KSTEP=3,  ALPHA(1)=0 .  DO,  0.D0,  -1.D0,  1.D0, 
T(l)=.  IDO,  TMAX=.  5D2, 

YZERO(l)=l . 0901751DO ,  ERR=.1D-11, 

H= . ID-1 ,  HMAX=  . IDO, 

1GFN=1,  IPC=1,  METHOD  =  l ,  INDEX*  0,  IAT=1 
$END 


User-Supplied  Subroutines 


SUBROUTINE  GFN  (G,H,F,Y,  J,T,  A,NR,NC) 
IMPLICIT  REAL*8  (A-H,  O-Z) 

DIMENSION  A  (NR,  NC),  Y  (4,  NC),  G  (NC),  T  (NC) 
G  (1)  =  T(J)  +  (l.-T  (J))*DEXP(-T(J)) 

RETURN 

END 

SUBROUTINE  AFNT  (A.N.T,  NR.NC) 

DOUBLE  PRECISION  A  (NR,  NC),  T 

A  (1,1)  =  -T 

RETURN 

END 
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An  Exact  Starter  -  BEGIN 


SUBROUTINE  BEGIN  (K,H,N,Y,  YN.YOLD, INVERT, IT,  METHOD,  IAT) 
PARAMETER  NR=20,  NC=20 
IMPLICIT  REAL*8  (A-H,  O-Z) 

EXTERNAL  INVERT 

COMMON  A  (NR,  NC),  ALPHA  (4),  T(NC) 

DIMENSION  Y  (4,  NC),  YN  (N),  YOLD  (4,  NC) 

DO  1  I  =  2,  4 
T  (I)  =  T  (1-1)  +  H 

1  Y  (I,  1)  =  DEXP  (-T  (I)  *  T  (I)/2.)  -DEXP  (-T  (I))  +  1.0 
RETURN 
END 


Source:  Krough  . 


Numerical  Results 


Give  results  at  t  =  1.0,  10,  and  30. 


t 

.10000000+01 

.10000000+02 

.30000000+02 


y  (t) 

. 12386512-01* 
. 12386512-01 t 

.99995460+00 

.99995460+00 

.10000000+01 

.10000000+01 


♦Produced  by  NLMS  methods  of  order  3  . 
t  Exact  solution 
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Problem  6 

y  =-100  t  y2;  y(l)  =  ^ 

Equivalent  Problem 

y'=  -100  y+100  y(l-ty);  y(l)=-^ 

Exact  Solution 

-  y(t)  =  -  o" 

1  +  f 


Classification 


(0,1,1,1» 


Eigenvalues  of  A 


{-100  \  . 


Methods  Applied 

NLMS  methods  of  order  3. 

Step  Size  Used 

Variable  step  size  with  initial  h  =  0.01. 
Time  Interval 

[1.  50]  . 

Special  Features  Exercised 
Variable  step  size 


0-  0  (t,  y). 


Exact  start. 
PCm  procedure. 
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FORTRAN  V  Inputs  Set-up 
Variable  step  size  ,  NLMS  order  of  3. 


$  IN  PUTS 

N=  1 ,  KSTEP=3 ,  ALPHA(1)=0. DO,  0.D0,  -1.D0,  1.D0, 
T(1)=1.DG,  TMAX=.  5D2 , 

YZERO(l)- .  19  6  0  7  843D-1 ,  ERR= . ID-1 1 , 

A  (1,1)  =  -100.  DO, 

H  =  .  ID-1 ,  UMAX  =  .IDO, 

IGFN  =  1,  IPC  =  1  METHOD  =  1,  INDEX  =  0 
$END 


User-Supplied  Subroutines 


SUBROUTINE  GFN  (G,H,N,Y,  J,  T,  A,  NR,  NC) 
IMPLICIT  REAL  *8  (A-H,  O-Z) 

DIMENSION  A  (NR,  NC),  Y  (4,  NC),  G  (NC),  T  (NC) 
G  (1)  =  100.  *Y  (J,  1)  *  (1.  -T  (J)  *  Y  (J,  1)) 
RETURN 
END 


An  Exact  Starter:  BEGIN 


SUBROUTINE  BEGIN  (K,  H,N,Y,  YN.YOLD, INVERT, IT,  METHOD, IAT) 
PARAMETER  NR-20,  NC«=20 
IMPLICIT  REAL*  1  (A-H,  O-Z) 

EXTERNAL  INVERT 

COMMON  A  (NR,  NC),  ALPHA  (4),  T  (NC) 

DIMENSION  Y  (4,  NC),  YN  (N),  YOLD  (4  (NC)) 

DO  1  I  *2,4 
T  (I)  -  T  (1-1)  +  H 
1  Y  (I,  1)  =»  1.  / (1  .+50.  *T  ffi) 

RETURN 

END 


Source :  Jain  . 
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Give  results  around  t  =  5,  10,  20,  30,  and  50. 


.5003  1250  +01 

.1000 1875  +02 

.2000  6875  +02 

.3000  4375  +02 

.5001  4375  +02 


.7983  6304-03* 
.7983  6304 -03 1 

.19988506-03 

.19988506-03 

.49963146-04 

.49963146-04 

.22215249-04 

.22215249-04 

.7995  3381-05 
.7995  3381-05 


Remarks 


The  decomposition  technique  is  applied  successfully  with  the  PCm 
procedure  for  this  nonlinear  problem.  Here,  9  (t,  y)  is  not  a  low-order 
polynomial,  but  the  program  can  be  seen  to  vary  the  step  size  automatically 
and  to  maintain  convergence  and  reasonable  accuracy. 

Problem  7 

y'=UZ  -  UBUy  ;  y(0)  =  (-1,  -1,  -1.  -1)  T 

where 


*Solution  produced  by  NLMS  methods 
t Exact  solution. 
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p 


W-  (Wj.  w2,  w3,  w4)  T  =Uy 
B  =  diag(  Pi,  02,  03,  04) 


=  1000, 

Exact  Solution 

y(t)  =  U 


p2  =  800,  fi3  =  -10,  P4 

f  01 

W-(l  +  Pj)  e  ^ 


1-  (1  +  P3)e 


0.001. 


02 

1-  (1  +  P2)  e 

04 

'V 

1-  (1  +  B4>  e 


Classification 


(0,  1,  1,  0). 

Eigenvalues  of  A 

{-1002,  -802,  8,  -2.001  }  initially 

{-1000,-800,-10,-0.001}  when  0.  001 1  >1. 

Method  Applied 

NLMS  methods  of  order  1 
Step  Size  Used 

0.01  <  t  <1,  h  =  .  ID-2, 

1  <  t  <10,  h  =  . ID-1 , 

10  <  t  <  1000,  h  =  .1D0. 


Time  Interval 


[o.oi,  100]  . 
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Special  Features 

Fixed  step  size. 

FORTRAN  V  Inputs  Set-Up 

To  maintain  better  precision,  matrix  A  (A(I,J))  and  vector  y 
(YZERO)  are  calculated  by  programs  specially  incorporated  in  MAIN. 

Fixed  step  size,  NLMS  of  order  1. 


$INPUTS 

N  =  4,  KSTEP  =  1,  ALPHA  (1)  =  -.IDO,  1.D0, 
T  (1)  =  .  ID-1 ,  TMAX  =  1.D3 
H  =  .  ID-2, 

IGFN  =  1,  IPC  =  0,  METHOD  =  1,  INDEX  =  0 
$END 


User-Supplied  Subroutines 


SUBROUTINE  GFN  (G,H,N,Y, J,T, A.NR.NC) 
IMPLICIT  REA*8  (A-H,  O-Z) 

DIMENSION  A  (NR,  NC),  Y  (4,  NC),  G  (NC),  T  (NC) 
DIMENSION  W  (4) 

WW  =  0.0 
DO  10  K  =  1 ,  N 
10  WW  -  WW  +  Y  (J,  K) 

DO  5  I  »  1,  N 
w  (I)  *  WW  -2.*Y  (J,I) 

5  W  (I)  s  W  (I)  *W  (I)  /4.0 
G  (1)  =  (-W(l)  +  W(2)  +  W(3)  +  W(4))/2.0 

G  (2)  «  (  W(l)  -  W(2)  +  W(3)  +  W(4))/2.0 

G  (3)  »  (  W(l)  +  W(2)  -  W(3)_+  W(4))/2.0 

G  (4)  =  (  W(l)  +  W(2)  +  W(3)  -  W(4))/2.0 

RETURN 
END 


7 

Source:  Krough,  Gear  . 
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Numerical  Results 


t 


Given  results  at  t  =  50,  500,  1000 


50 


-.50095236+01 

-.50095561+01 


-.500953  2  6+  01 
-.50095561  +  01 


.49904764+01 
.49904439  +01 


-.49904764+01* 

-.49904439+01t 


500 


-.50007681+01 

-.50007688+01 


-.50007681+01 

-.50007688+01 


.49992319+01 

.49992312+01 


-.49992319+01 

-.49992312+01 


1000 


-.50002903+01 

-.50002905+01 


-.50002903  +  01 
-.50002905+01 


.49997097-+01 

.49997095+01 


-.49997097+01 

-.49997095+01 


Remarks 


Since  only  three  different  step  sizes  are  used,  there  are  only  three 
matrix  inversions  and  three  matrix  exponentials  needed.  To  reach  t  =  1000 
takes  about  4  min  46  sec  with  a  maximum  error  ~  .38  x  10-7. 


♦Solution  produced  by  NLMS  methods. 
tExact  solution. 
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9.  CONCLUSIONS 


A  family  of  strongly  stable  and  consistent  NLMS  methods  has  been 
developed.  When  applying  NLMS  methods  to  the  problem  y  =  Ay  ,  which 


implies  g(t,  y)  =  0,  NLMS  method?  produce  the  approximate  solution 

y^  =  e  yo  =  e  V() .  Since  the  matrix  exponential  is  computed  by 

the  Pade  approximation,1**’ 1  which  is  stable,  the  lim  ||  y  ||  — >  0  as  n  a> 

13  n 

establishing  the  A-stability  in  the  sense  of  Dahlquist. 


If  9  (t,  y )  is  a  slowly  varying  function  that  can  be  approximated  by  a 
low-order  polynomial  in  t  and  A  is  a  slowly  varying  function  in  t,  NLMS  methods 
allow  the  use  of  a  larger  step  size  and  can  produce  almost  exact  solution  in  the 
absence  of  ruund-off  errors.  In  the  event  that  g  (t,  y  )  is  not  a  low-order 
polynomial  in  t  or  not  a  slowly  varying  function  in  t,  NLMS  methods  can  still 
produce  accurate  results  if  not  too  large  a  step  size  is  used. 


In  general,  because  of  the  initial  local  discretization  error,  implicit 
NLMS  methods  are  better  than  explict  NLMS  methods.  High-order  NLMS 
methods  are  better  than  low -order  NLMS  methods.  In  the  event  that 

g(t,  y )  is  a  constant  in  t  and  y  ,  we  see  that 


1=0  i=0 


Therefore,  explicit  NLMS  of  order  K  is  equivalent  to  implicit  NLMS  of 
order  K.  In  this  case,  low-order  NLMS  methods  will  suffice;  high-order  NLMS 
methods  are  not  necessary.  In  fact,  high-order  NLMS  methods  require  more 
computations,  a  good  Pade  approximation,  and  an  accurate  matrix  inversion. 

Several  numerical  experiments  has  been  performed  upon  test  prob¬ 
lems  utilizing  different  selection  of  a  .,  from  both  strongly  and  weakly 

stable  NLMS  families.  The  numerical  results  agree  to  nine  significant 
digits. 
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The  numerical  results  presented  in  this  report  mainly  are  used  to 
demonstrate  that  the  NLMS  program  can  produce  satisfactory  results.  The 
computations  herein  were  intended  for  illustrative  purposes  only.  It  is  our 
intention  to  apply  NLMS  methods  to  more  difficult  problems.  We  intend  to 
develop  criteria  to  determine  how  good  NLMS  methods  are.  These  will  be 
reported  separately. 


The  NLMS  package  requires  the  user  to  guess  an  initial  step  size  h. 
However,  the  package  is  "robust";  it  determines  the  most  favorable  step  size 
for  rapid  computation.  Experience  indicates  that  less  total  computation  time 
results  if  h  is  chosen  too  small,  rather  than  too  large. 

If  a  large,  stiff  system  is  to  be  solved  with  the  variable-step-size 
technique,  an  eigenvalue-eigenvector  technique  such  as  EBSPACK  is  recommend¬ 
ed  to  replace  PADE.  This  can  save  a  large  amount  of  computation  time. 


58 


TR  5341 


10.  REFERENCES 

1.  D.  Lee,  Nonlinear  Multistep  Methods  for  Solving  Initial  Value  Problems 
in  Ordinary  Differential  Equations,  Ph.  D. dissertation,  Polytechnic 
Institute  of  New  York,  1974.  Also  NUSC  Technical  Report  4775,  24  May 
1974. 

2.  P.  Henrici,  Discrete  Variable  Methods  in  Ordinary  Differential 
Equations,  John  Wiley  and  Sons,  Inc. ,  New  York,  1962. 


3.  H.B.  Keller,  Numerical  Methods  for  Two- Point  Boundary -Value 
Problems,  Blaisdell  Publishing  Co.  ,  Waltham,  Massachusetts,  1968. 

4.  J.  Certaine,  "The  Solution  of  Ordinary  Differential  Equations  With 
Large  Time  Constants,  "  in  Mathematical  Methods  for  Digital  Computers, 
A.  Ralston  and  H.  Wilf,  eds. ,  John  Wiley  and  Sons,  Inc.  ,  New  York, 
1960,  pp. 128-132. 

5.  L.F.  Shampine  and  M.  K.  Gordon,  Computer  Solution  of  Ordinary  Differ¬ 
ential  Equations.  W.  H.  Freeman  and  Co.  ,  San  Francisco,  1975. 

6.  B.L.  Ehle,  A  Comparison  of  Numerical  Methods  for  Solving  Certain 
Stiff  Ordinary  Differential  Equations.  Department  of  Mathematics, 
University  of  Victoria,  Canada,  Report  70,  1972. 

7.  C.W.  Gear,  Numerical  Initial  Value  Problems  in  Ordinary  Differential 
Equations.  Prentice -Hall,  Inc.,  Englewood  Cliffs.  New  Jersey,  1971. 

8.  R.K.  Jain,  "Some  A-Stable  Methods  for  Stiff  Ordinary  Differential 
Equations, "  Mathematics  of  Computation,  vol.  26,  no.  117,  pp.  71-77, 
1972. 

9.  F.T.  Krough,  "On  Testing  a  Subroutine  for  the  Numerical  Integration 
of  Ordinary  Differential  Equations,  "  Journal  of  the  ACM,  vol.  20, 
pp.  545-562,  1973. 

10.  L.E. Weaver,  Systems  Analysis  of  Nuclear  Reactor  Dynamics,  Rowman 
and  Littlefield,  Inc.,  New  York,  1963. 


59 


TR  5341 


11.  R.A.  Willoughby,  Stiff  Differential  Systems.  Plenum  Press,  New  York 
and  London,  1974. 

12.  J.  L.  Blue,  and  H.  K,  Gummel,  "Rational  Approximations  to  Matrix 
Exponential  for  Systems  of  Stiff  Equations,"  Journal  of  Computational 
Physics,  vol.  5,  1970,  pp.  70  -83. 

13.  G.  Dahlquist,  "A  Special  Stability  Problem  for  Linear  Multistep 
Methods,  "  BIT,  vol.  3,  1963,  pp.  27-43. 


60 


Appendix 
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the  characteristic  coefficients,  alpha 
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METHODrO  —  LMS»  METH0D=1  -*  NLM5 
FiXED-STEP-bIZE  IV50,  INDEX=0  PREDICTOR 

IVxO j  index=i  corrector 

VARIABLE-STEP-SIZE  IV  NE  0#  PREDICTOR-CORRECTOR 


58  IF( METHOD, NE.O)  60  TO  64 

59  IF(INOEX,EQ.O)  CALL  LMS<KSTEP,H.YOLD,N»YN,IZ£RO) 

60  IF  (  IV,£Q,0,  AND;  INDEX  aEQ*i  )  CALL  LMS  (KSTEP,H,  YOLL)  »N,  YN»  I  ONE ) 

bl  60  TO  59 


TR  5341 


* 

* 

* 

« 

* 

* 


*  *  *  # 
*  » 

*  « 

*  * 

H 

*  * 

4 

*  * 

i-i 

o 

* 

*  * 

«• 

oc 

Id 

#  « 

o 

Id 

2 

*  * 

H 

N 

O 

H 

»  * 

•k 

*4 

i-H 

4 

*  # 

ac 

m 

•4 

*  # 

cc 

z 

z 

m 

«  • 

o 

H 

>- 

13 

*  # 

»■* 

o 

m, 

* 

M 

*  « 

Id 

H 

z 

2 

*  * 

_l 

* 

«> 

Id 

*  * 

9 

^  UJ 

Cl 

O 

Z 

*  » 

H 

UJ  z 

_l 

-1 

o 

*  * 

9 

z  o 

o 

«o 

n  H 

#  * 

9 

o  *-• 

>-  >-  >-  1-  » 

j>  *• 

o 

•-I  » 

«b 

<  - 

<  o 

*  * 

0. 

•  at 

I  M  I  w(t 

«  * 

*  UJ 

<4 

»  44 

•  Ui 

*•  * 

UJ 

UJ  z 

Q.OQ.ON 

«  « 

u 

Z  H 

UI  H  Ui.H  *4 

jt  jk 

z 

H  » 

H 

•  H 

4b  ^ 

4*  * 

UJ 

•  z 

IAZUIIZ 

*•  #■ 

9 

z  • 

5Z  M  *  t-t  >- 

*  CM  » 

DC 

•  o 

*9  « 

UJ 

Q  _J 

9 

l/) 

Z 

H-UJU.  * 

> 

J  O 

x 

X 

• 

»x  o  *- 

z 

O  H 

-J 

_l 

Q  O' 

•#  H  » 

o 

H  • 

2 

z 

-1  <o 

*  HOC  * 

u 

•  I 

a 

#•  o  .* 

I  • 

-1 

-1 

>-  o 

#  I0H  * 

o 

•  0- 

_i 

J 

•  h 

*  u  # 

~z 

(L  UJ 

< 

< 

X 

#  H  4  * 

o 

UJ  H 

o 

u 

•  o 

#Ult  # 

H  • 

h  in 

CL  9 

Hid  *• 

9  H  rlZ 

9  2 

«"* 

UJ 

*  K  4  * 

H  •**»• 

o 

o 

H  — 

jtae  .# 

Z  -X 

-*in 

• 

• 

1/)  O 

•  oh* 

O  *-•  0-  X 

9  X 

a 

ui 

X  •O' 

•  urn* 

H  X  — 

IJO 

UJ 

z 

—  UJ® 

«  « 

Z  Q 

JZ  H 

* 

• 

9  -  z 

*  l-Q  # 

O  MIIUJUJ 

X 

K 

X  M  »D 

#  9  J  # 

9  X  Z  H  X  a. 

J  JO 

Id 

hi 

-1  ^  H  H 

toe  *• 

XHMUO 

J  _J  H 

o 

I 

z  z< 

•  X  Z  • 

«•  it  x  at  <  h 

<  4 

2 

KH.a 

*  4.  « 

H  XH  Old  9 

u  u  o 

►, 

HO 

-1  II  9 

»HI  # 

X  X 

9 

•s 

J  **  ► 

#  <  U  # 

J  ^  ^  ^ 

Q 

o 

<  hq:  ~ 

*  « 

ZZIOXIO 

HI  O  O 

z 

z 

CJ  •  O  *4 

*  QCUJ  • 

H  •  HH  X 

•  •  H 

< 

< 

H  • 

H 

•  ON* 

•  UJ  X  x  «o>  • 

O  O  K 

• 

• 

«z  ♦  *z 

•  H  M  * 

DC  _l  X  X  H  >0  H 

Z  DC  UJ  O  • 

o 

o 

o  •  a.  o  9  • 

z 

•  CM/)  • 

ac  •  9-90 

%  m  •  •  m 

• 

• 

•  d  UJ  •  ji  H  V 

»  d  1  * 

o  •  •  »X  *ZHN  QQZ 

9 

0 

UJ  II  H  Ul  •  II 

ii  ai  •  ac  o.  » 

OK  O  HH  IO  UJ 

IIHOOO  • 

Id 

UJ 

Zh^Z  h 

4>  •  oc  uj  tOHa:ajJZ''ZNHiiioxx 

• 

• 

•  SC  »H 

ho 

•  OH# 

It  IIO  •  •  .HHHH 

«  H  H  H 

> 

> 

>  ®  ~»X  ® 

o  •ui/)«a:xuAi  x<xii*4»-«ouju)x 

H 

M 

M  ®  Q  M  _J  ® 

* 

h*  *  a  x  *-•  v  x  x  +4  x  >4 

H  >HXX  J 

«■> 

w 

<«  |  w  w 

uJ 

H  ««»  <• 

u. 

u. 

U.OOU.U.OZO*  * 

UOU.II  U.  Ul  U-  O  U.  w 

0"0lllUl 

*4 

»-4 

H  3  >-HHQ 

>- 

9  *  *  #  ♦ 

HH  HI  HHHIL 

*  o*  > o  «  O'  o  o  *4  m 

•o  \n  \ o  so  ®  h  h  r»  h 


o  uuu 


tMi00,9®H-®0‘O*4AI»03,in<0N®0'O*-IC¥»00,9®H®9OH(\|i0 

■o®®®®®®®HHHHr>*HHHHH®®»«®®®®®®9O'0'a* 


A-35 


'mmmm 


tarmi 


T R  5341 


z  - 
*z 

H  >- 

II  H 

K 
^  *- 
Is-  X 
_  Ui 
02 
Q  *- 

N 

N 


*  *  « 
*  * 
*  * 
*  # 
*  » 

;  * 
*  * 

*  # 
*  * 
*  # 
*  * 

*  « 

•  » 

*  « 

*  * 

*  * 

*  # 

*  « 

*  *  * 

*  £C  « 

*  O 

*2  * 

*  # 
♦•x  *. 

#  <  # 

»#  r* 

*>  «. 

*  ©  # 

*  2  # 

I  N  | 

#  W  # 

*  # 

*  *  * 

*  Ui  # 

*  o  » 
*2? 

*  Ui  * 

*  o  • 

#-cc  # 

*UJ  # 

■# 

*  2# 

JO  # 

*  # 

*u»  # 

#  or  # 

#  o» 

#  »-  # 

*  o  # 

*  UI  # 

*5  # 

*  QC  « 

_  *  O  # 

it 

sBBiSi 
*  * 
ooz  2  • 

<°u<2*  2 


+ 

a 

Ui 

i- 

in 

x 

w 

a 

Jo 
O  00 

X 

v  o 


">o 

*  ® 
UI  ^ 
2  ~ 
>-  -> 

•  ■) 

•>  © 

«k  W 

?S 
&2  = 
t:  .2 

,  *  ►-  a 

^ O  w 
*°  •  L0 


» 

# 

* 

# 

# 

* 

* 


* 

#■  * 

*■  ^ 


*## 
* 
* 
« 
* 
* 
# 
* 
* 
* 
* 

* 

# 

* 

# 

* 

# 

* 

Jf 

^  ; 


<>i 

® 


«  X 

$  ir 

®  -N- 


© 

K 


«£»  3 

*  >“  QC  Q 
*-  O  II 
®»2J 
®  ^  <  CC 

a  u>  5$ 


ar  >h 
£?£ 

f  §■>«“ 

"j»:0 

J:«3° 

P  U.  O  O  o 

u  M  J  >-  to 


o 

« 


o  o  t> 


<v  -» 
®  ® 


* 
# 
» 
*■ 
«• 
* 
* 
# 
# 
* 
» 

« 

* 

.# 

« 

# 

* 

-*  # 

•"*  # 

ir * 

*  # 

UJ  « 

5  * 

>-  * 

II  ft 

2  ft 
*  *# 

II  ft  * 

* 

Ui  # 
nt.  I 
®in  # 

_  x  * 

O  w  ji 

o>-  # 


OiiO 
®  ® 


■f  r 
a.  j? 
uiin’t 

J* 
o« 

i-in  * 
>*  ui  ♦ 

9l  # 

**  2  * 

M  * 

~Sri 
*°# 

*“.*  ft 

2^: 

~  ^  * 

» 

££  # 

■K  O  # 

Wn» 

w._  # 

^2: 

^  OL  # 

5?  -i  * 
a:  u  # 

» 

JL  * 

#  #  * 


a; 

ui 


2 
§ 

*  k 

Ui  UJ 
?  K 


Ui 


■5 

>» 

II 


o 

11 

»■*  X 
II  Ui 
h-  Q 
2  2 


(9  a 
UI  uJ 


X 

<  a. 

2  ui 

»“  i- 

in 

•S22 


M  II 

1  H 


2  (o 

Hkt  - 

ae  o 
Q-  2  < 
.^•U itnin 

2  2  **■ 0  o 

U  u  H  Q  3 


♦ 

® 


o 

X 
M  UJ 
*  ISI 
-<  V 

■f 

■>  ~ 

^  »■* 

>-  •  Ui 

II  a  o 
-*-UJ  z 

i2* 

*•  ~*8i 

in 

® 


<v 


?Ooo 

•  •  a, 
r  <3  O 
{7  UJ  UJ 

•io^t 

-*  «  M  J 

^  ui  —  3d 

a.ii. 


L> 


Uy  y 


-'=ss^3332333S5533335Sa! 


T 


< 

I 

! 

I 

J 


I 


TR  5341 


ft  *  * 
*  * 

ft  O  ft 
ft  Z  # 
ft  *4  * 
ft  3  ft 
ft  3  * 
ft  3  ft 

*  O  « 

ft  Q  * 
ft  ft 

ft  13  ft 

»  r  ft 
ft  o  * 

ft  JL  * 

*  UJ  * 

*  od  * 

*  * 
*  Ji  * 
ft  U  ft 
Jt  Z  .* 

ft  -«  ft 
ft  K  ft 
ft*  ft 

ft  fO  ft 
■*  .* 

ft  oc  *- 
ft  O  ft 
ft  li.  ft 
ft  ft 

ft  >-  ft 


ui  *- 
»-  II 

II  o 
-*  ac 
h  ui 
*-N 


a. 

UJ 

►- 

i/> 

* 

2  — 
*  V- 
H  II 

’  II  —* 


<o  o 

©  X 
ul 
O  Ni 
3  >- 


« 
■ft 
ft 
ft 
ft 
ft 
ft 
*■ 
ft 
.* 
ft 
ft 
ft 
ft 
ft 

—  ft 
►-»  ft 

*  ft 
»-«  ft 
ft  ft 

a. » 
ui  » 
k  * 
co  » 
*  * 
«  * 
>-  » 
ii  « 
» 

»-•  * 
»  ft 
*4  ft 

—  ft 
>-  ft 


H-  ft 
ZJk 

<  ft- 
»-  ft 
l/l  ft 
Z  ft 
O  ft 
U  ft 

ft 

x*» 

ft 

3  ft 
3  ft 
ft 
ft 
IA  ft 
UJ  ft 
3  ft 
3  ft 
3  ft 
in  * 

ft  -4 

z  *  ♦ 

►4  ft  Q. 

<  ft  Ui 
I-  ft  b~ 
z  *  to 

—i  +  —t 

<  *  ii 
Z  ft  Q. 

ft  Ui 
ft  H- 
•  UI 
ft  ft  *■« 


< 

X 

X 


O  ft  O  3 
»-  II  o  • 

i/»  a.  a  o 

HU  .Ih 
H  H  i\|  -  Q 

li.lAIIU.OZ 

HHtHlIU 


CL¬ 

UJ 


* 

* 

* 

ft  ft  ft  ft 

ft  ft 

ft 

M 

ft 

ft 

c_> 

c  z 

ft 

ft 

CM 

-  >~ 

* 

ft 

CM  - 

ft 

ft 

o 

I  >- 

ft 

ft 

iM 

«. 

ft 

ft 

w 

x  3 

* 

ft 

»- 

-  X 

ft 

ft 

o  *• 

ft 

ft 

3  — . 

-  m 

ft 

ft 

3  O 

X  3C 

ft 

ft 

-  M 

<JC  *■ 

ft 

ft 

— .  * 

lO  CM 

ft 

ft 

CD  O 

U  ,lc 

ft 

ft 

>M  M 

*»  ^ 

ft 

ft 

o 

m  « 

X  -4 

ft 

ft 

C\J 

O  3 

4  i 

ft 

ft 

«•  o 

CM  *-* 

CM  - 

ft 

ft 

O  C\J 

<*•  * 

UJ  ft- 

— * 

ft 

ft 

CM  — 

«■  M 

K 

ft 

ft 

*■»  e> 

>— <*  » 

T-  r* 

< 

ft 

ft 

3  * 

M  — > 

•*  3 

H 

ft 

< 

ft 

T  «■* 

X  o 

UJ  - 

ft 

X 

ft 

<  o 

0  CM 

^  3 

1- 

ft 

10 

ft 

»  CM 

h-  1- 

H 

ft 

o 

ft 

— % 

Cl 

4  - 

ft 

X 

ft 

o  o 

o  *M 

rO 

in 

ft 

a. 

ft 

CM  CM  CM  * 

<  K- 

H 

ft 

in 

ft 

*  H-r 

*  ct 

X  - 

«u 

ft 

2el 

-H* 

HI 

ft 

o  X 

o  ^ 

a.  cm 

X 

ft 

-J 

•— 

ft 

CM  < 

CM  fO 

3  i- 

UJ 

ft 

h 

<t  u 

2 T 

ft 

*~  K) 

*  < 

<  * 

ij 

ft 

Or 

X  z 

>- 

ft 

IO  LJ  J  « 

»  ■— 4 

ft 

2: 

< 

ft 

X  * 

—  O  J  h 

M 

ft  X 

1- 

z 

ft 

<  <* 

MUM 

r  - 

ft  h- 

1/)  |-<  D 

H 

ft 

»  O  i  (\l  - 

<3  4- 

ft 

M 

o 

ft  — 

<-»  CM 

x  *  r 

•b  » 

>-  *  x  x  >  x  x 

**OOUD~ 
Z  ft  O  Z  3  *4 


0*00  *  -  O  >-  r<->  >-• 

ujfttMtMo^fM.rr 

X  ft  *-  »tVI  Ci  <*  ~  C  0 


CM 

ft 

jo  Hi/ia 

o 

ft 

1-0*' 

CM  3 

o 

«ft  «b 

>- 

ft 

4UX 

K 

ft 

rvi  X 

ft— 

CM 

"M  «-< 

ft 

IL  U  O  • 

—ft 

3  l/l 

ft 

-» 

w  < 

O  >M 

«b 

X  X 

o 

W 

X 

#1H3  • 

U  Z 

ft 

* 

CM  CM  CM 

s 

<  X 

1/ 

& 

* 

ft 

UJ  3  i/1  • 

z 

*  in 

ft 

w 

X  ui 

ft^> 

<— 

«  «b 

r™ 

«* 

cx 

ft 

1-  >-  J  » 

•  >-« 

ft 

<t 

H 

—ft 

X  •-* 

ft 

JL 

UJ 

ft 

i/1  >-  3  3  I 

• 

» 

ft 

X 

»■  —ft 

a. 

o 

<t  a. 

w 

o 

JL 

H 

►- 

ft 

m.  OD  <t  <  < 

• 

•X 

ft 

X 

-ft  o 

CM 

—ft 

C/1 

ft 

»-  UJ  3  — 

• 

-  o 

ft 

3  O  CM 

—ft 

«b 

O  <t  (X 

Am 

< 

ft 

J  OC  1- 

<t  I- 

* 

■"M  * 

o 

M 

■A 

ft 

ft 

3  UJ  J  Z  X  0 

*^3 

ft 

«b 

*  O  !\l 

CM 

•k 

W 

10 

ft 

5.  3  <  <  uJ 

ft-ft 

UJ  uJ 

ft 

—ft 

O  CM 

o 

o  o 

5 

#  3 

X 

>  z 

o  > 

ft 

C>  CM  —  O  ~t 

CM  *h  -« 

3 

*  X  < 

i/1  0 

z  u. 

< 

ft 

-M  *-  r  CM  —> 

—  l/l 

Z 

ft  c  u 

H-  O  i-»  3  X 

A 

ft 

%  X4  ^  H 

#  3 

Z  X 

o 

ft 

O  <  UJ  X  ■* 

<30 

UJ 

*  z 

UJ  X 

M 

ft 

:M 

3  '3 

-w 

ft  ‘H 

z 

t- 

ft 

*-  z  z  »  z 

*  <r.  t 

#  3 

3  J) 

3 

ft 

<  o  o  o  o 

O  X  X 

»- 

#  z 

0  l-l 

3 

ft 

H  H  W  H 

H 

3 

ft  3 

■X  X 

O 

ft 

Z  CO  iO  J\  A 

ji  j  j 

O 

#  X 

<r  k 

CO 

ft 

a  ^  2  z  ? 

Z  J  3 

X 

ft 

ft 

X  UJ  3  UJ  3 

3  3  3 

3 

ft 

ft 

Z  -£  X  X  € 

X  3  3 

3 

ft 

ft 

3  M  M  M  -1 

-«  3  3 

i/1 

ft  ft  ft 

«  * 

ft  ft 

ft 

ft 

ft 

3  3  3X3 

3X3 

<0 

3) 


U  JU 


O JOUUUJUUU 


-4fMO^ioor*-coo‘o*Hoji04>in 

(OiO'O'OiOiOO'OiOt 


H  H  H  H  H  H  H  H  H  H 


A-37 


I 


TR  5341 


« 

*  *  *  * 

ft 

• 

ft 

* 

ft 

• 

ft 

* 

z 

ft 

ft 

u. 

ft 

* 

ft 

» 

H 

ft 

* 

X 

ft 

ft 

a 

ft 

* 

ft 

* 

Li 

ft 

* 

C/1 

ft 

* 

or 

ft 

» 

UI 

ft 

• 

> 

ft 

* 

z 

ft 

• 

H 

ft 

* 

ft 

* 

< 

ft 

* 

ft 

* 

fe. 

ft 

* 

H 

ft 

• 

x  a 

ft 

* 

Ul<  u 

ft 

* 

m 

ft 

« 

X  C/l 

ft  X 

* 

in  x  < 

ft  Ui 

ft 

HliJCL 

ft  H- 

* 

>- 

*  in 

« 

z  -  u 

ft  it 

ft 

(il  X 

ft 

* 

I  <  Ul 

ft  - 

« 

3t  X 

ft  -- 

ft 

in  < 

ft  o 

• 

Z  Ui 

ft  o 

ft 

o  >  in 

ft  lO 

ft 

ft  - 

« 

Kino 

ft  o 

# 

<  M 

ft  o 

• 

m  a  »- 

ft  <M 

* 

*-•  z  <t 

ft  - 

ft 

J  4  3 

ft  o 

« 

4  3 

ft  o 

* 

a:  ►*  in  i_> 

ft  H 

* 

OHUJ  J 

ft  — 

#  *-•"!  1-  4 

ft 

*<2<U 

ft  o 

ft 

OM  J 

ft  1- 

*  -•  3UI 

ft 

*  3  i/»  O  > 

ft  o 

«ZUI  JO 

ft  o 

* 

►ft  O  4  03 

ft 

» 

0  0  4 

ft  »■* 

*  z 

ft  -« 

ft  4  X  X  *-« 

ft 

» 

4  4  • 

ft  • 

#  jut  tH 

ft  »- 

*►-•003 

ft  Id 

* 

o  o  • 

ft  • 

tvflir  ii/» 

ft 

*  •-«  x  x  ►* 

*  in 

» 

*  ►■« 

« 

*  — 

* 

* 

#  ft  *  * 

ft  X 
ft  ►ft 

o 

a 

z  - 
•  o 
il 

il  — 
3  3 


O  rH 

j  a. 


O  O  O 
O 033000 

• Q  •  »Q  o  j 

o  •  o  o  •  •  • 

II  O  ||  IIOOOO 

—  II  «-  *■»  II  II  II  Q 

3  3  3  O 

-  3  -  *333  • 

H  »  H  H  »  <k  k  H 
^  H  -  -  MH  M  II  ^  2 

"IX - -  —  ►  «• 

Z  4  <M  *1  I  X  I  -II  II 


o 

3 

• 

O 

II 

z  — 

»  it 

H  k 

II  “> 
at  - 


o 

Q 

• 

O 

II 

—  z 
*  - 
— » 
II 


030303 


O 

Q 


»-  IQ  <»  IO  —  —  «o 
»-t  »-«  >-* 

Z  O  O  O  X  X  o 
3  3  3  3  0.  30 

CM  -ft  K1 


X 

LJ 

t- 

J) 

it 


o 

O 

*■*  «o 

3  - 
-o 

H  O 
(\J 

4  - 

Z  *  o 
-  X  o 

k-H  II  — *  I 

II  ~~ 
3  3 

-  O  ■ 
<0  i 

0X0 
3  4  3 

3 


o  * 


o 

u 


ft  ft  ft 

« 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft  z 

ft 

ft  o 

ft 

ft  ™ 

ft 

ft  1- 

ft 

*  4 

ft 

ft  •  M 

ft 

ft  Ui  ►■• 

ft 

*  a:  _i 

ft 

ft  Ui  < 

ft 

*  x  >-• 

ft 

*  »- 

ft 

*  in  ►•« 

ft 

*  ►-  z 

ft 

ft  X  >-* 

ft 

ft  4 

ft 

*  i—  in 

ft 

*  in  ui 

ft 

*  o 

ft 

ft  X  a 

ft 

ft  Ui 

ft 

*  t-  z 

ft 

*  in  o 

ft 

ft 

#  i-  i- 

ft 

*  JO 

ft  Z  O 

ft  3U1 

ft  -  3 

*  x  in 

ft  H  • 

o 

• 

ft  II  o 

o 

ft  X  id 

ft  3  11 

• 

ft  4  Z 

ft  Z  ~ 

o 

ft  UJ  ►* 

ft  -  KJ  3 

II 

ft  z  z 

*  IO  - 

o  — 

ft  M  z 

ft  II  H  -•  3  H 

ft  J  ►-• 

ft  -»  — 

•  «h 

ft  Z  id 

ft  O  kft 

O  H 

ft  O  Ui 

ft  M  i  J  X 

II  - 

*  z  a 

ft  to 

~  OJ 

ft  -* 

H  w 

• 

ft 

H 

* 

ft  o 

j:  r 

«  »  « 

ft  3 

m 

r) 

—4 

>-  X 

0  0  3  3 

OH^l04’irt^M0JOH^i0iinfiM00'Ort[y|i0^lfl>0h#3‘OH.M^*in4) 


A-38 


i 


I 


-  * '  1 


TR  5341 


* 

ft 

ft 

ft 

ft 

ft  ft 

ft 

* 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

* 

ft 

ft 

ft 

* 

ft 

ft 

ft 

* 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

0"* 

ft 

ft 

ft 

ft 

ft 

ft 

« 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

* 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

-> 

ft 

I 

ft 

ft 

ft 

•> 

ft 

< 

ft 

ft 

ft 

H 

ft 

- _ 

ft 

ft 

ft 

ft 

a 

ft 

ft 

ft 

H 

ft 

X 

ft 

ft 

ft 

%— 

ft 

UJ 

ft 

ft 

ft 

ft 

ft 

ft 

— » 

ft 

X 

ft 

l/l 

ft 

ft 

-ft 

ft 

CL 

ft 

2 

ft 

ft 

ft 

ft 

* 

-ft  O 

ft 

M 

ft 

ft 

2 

ft 

I 

ro  r- 

ft 

< 

ft 

ft 

«— 

ft 

* 

-ft  -ft 

ft 

y- 

ft 

ft 

>- 

ft 

3 

ft 

2 

ft 

ft 

-* 

ft 

% 

o  o 

ft 

O 

ft 

ft 

o  _J 

ft 

u. 

H  1- 

ft 

o 

ft 

ft 

►  4 

ft 

V— 

•* 

ft 

ft 

ft 

-ft  ~z 

ft 

X 

w 

o  o 

ft 

X 

ft 

>-* 

#  Ui 

—  M 

ft 

4 

>- 

O  3 

ft 

4 

ft 

% 

H 

ft 

O  *-•  X 

ft 

ft 

* 

ft 

Uj 

ft 

o 

H 

ft 

O  X 

ft 

— ► 

A  — » 

ft 

ft 

z 

ft 

X  X  10 

ft 

o  -ft 

ft 

* 

ft 

I 

D 

ft  1- 

UJ 

ft 

•k 

•  • 

ft 

A 

ft 

< 

♦ 

ft 

UI  l/l  J- 

ft  o 

a 

H 

T»  C5 

ft 

-ft 

ft 

m 

— % 

ft 

T-  UJ  3 

•  in 

w 

U  ui 

ft 

1 

ft 

ui 

■>  * 

1-  X 

ft  -ft 

-ft  ~~ 

X 

•  • 

ft 

ft 

^  *> 

ft 

X  4  X 

ft 

O' 

X  - 

< 

X  X 

ft 

ft 

ft 

X  X 

H 

ft 

UJ  -J  o 

ft  o 

o 

1  X 

UJ 

jj  jj 

* 

ft 

ft 

ft  4 

* 

1-3  0 

ft  1- 

-ft 

—  jj 

♦ 

O  3 

ft 

ft 

O  -ft 

~  IO  ~  CM 

l 

ft 

l/l  u 

# 

~)  »- 

2  2 

ft 

X 

ft 

CM  - 

-ft  c\|  -4  UJ 

4 

ft 

1 

_l  o 

ft  o 

o 

►  Ui 

H 

►ft  M 

ft 

< 

ft 

H  H 

CL  -ft  *  » 

UJ 

ft 

-ft  4  -ft 

ft  o 

*- 

►ft  ftC 

ft 

w 

ft 

•> 

»  •  H  X 

ft 

ft 

1 

O  -ft 

ft 

■»  — 

X 

•  • 

ft 

ft 

O  X 

2  lO  —  < 

— * 

ft 

X 

ft  — 

o 

-ft  >- 

r 

OQu 

ft 

l/l 

ft 

CM  <t 

►  CM  4  UJ 

-ft 

ft 

-j  i noc 

ft  C  J  U> 

— •  — i 

CL 

0  2  2: 

ft 

2 

ft 

H  N 

X  -ft  — 

■W 

ft 

2 

OO 

ft 

►ft  2 

-J 

•  <  4 

ft 

-ft 

ft 

-ft  3 

4  -ft  X  c 

4 

ft 

-ft 

ft  • 

x  ► 

< 

o  •  • 

ft 

< 

ft 

*  a 

—  *  X  •/  ’ 

X 

ft 

1- 

CO 

ft  UJ 

O  Z 

T 

X  X  2  2 

1 

2! 

II 

ft 

b- 

ft 

o  o 

»-  lO  UJ  <  -  - 

X 

ft 

►ft  1  o 

ft  2 

•ft 

* 

II  ►  ►  - 

* 

~  H  -1 

ft 

2 

ft 

CM  • 

ft  CM  J  —  -ft  -ft  J  ft 

u  o  — * 

ft  • 

•  H 

ft* 

ft* 

O  -ft  -ft 

-* 

•4  •  • 

ft 

O 

ft 

H  -ft 

cm  uj  -ft  ii  in  uj  ii  ii  < 

ft 

H  o 

ft 

r*  II 

II 

II  3  —  II  II 

w 

M 

II 

►  1—  1— 

ft 

O 

ft 

II 

<M 

>  —  cm  3  «-•  “5 

11 

ft 

J  J  X 

ft  X 

UJ  -ft  3  'SC 

►  2  *ft  O 

z 

H 

m'}8 

ft 

ft 

-ft 

X  —  -ft  -ft  4 

ft 

X 

o 

•  UJ 

• 

►ft  X 

►  •  • 

ft 

H 

ft 

-ft  H 

*— •  -ft  •  X  "1  Cl 

ft 

x  o  o 

ft  Q 

lO 

n  in 

►  30-0 

II 

z 

o 

o 

K)  |/|  ^ 

ft 

CL 

ft 

1  * 

o 

1  -ft  O  o  O 

ft 

UJ  3  -J 

ft  2 

►-  O  O  O  -ft  O  o 

* 

«  W  H 

ft 

ft 

2  -# 

-J  2  —  J-  _J  -ft  -ft 

.* 

ft 

ft  -ft 

M  H 

-ft 

-ft 

-  JHH 

»* 

J  H  H 

H  —  — • 

ft 

ft 

_J  —  X  J 

ft 

ft  — 

►ft  J 

w 

H-* 

I  UL  U. 

ft 

ft 

X  H 

o 

4u.40400x 

ft 

ft  X 

u.  o 

0  0X403 

z 

J 

o 

o 

Q.  -ft  -ft 

ft 

ft 

ft 

►ft  X 

o 

J  -ft  aJ  5  'J  J  3 

1 

ft 

ft 

ft  ft 

ft  -ft 

►ftooaxoo 

3 

CM 

O 

O  CM  "O  10  Si 

IO 

in  -ft 

CO 

O' 

>o 

M 

CM  M  CM  >vi  M 

o 

O  IO 

o 

c 

-ft 

-ft 

—ft  -ft  -ft  —ft  *ft 

.ft 

-ft  -ft 

* 

«* 

u  JU 


U  JUUU 


M8}>0  H  M  O  J- 

intnn.o-o-o.o.oooo.ooi-f^rT-r'-f-r-r-i^r'-.oo.o.o.o.oo.o.o.oJMJ'J'j* 


A-39 


1 10  YU  ( 1 )  -  YN  ( 1 )  -mL^A  ( 1 )  *EAH  ( I  *  J )  ♦  i  <  1 ,  J ) 
RETURN 


TR  5341 

r 


\ 

| 

H 


*  *  *  * 
* 

# 

* 

« 

« 

* 

* 

« 

* 

* 

* 

* 

*  < 

*  v 

* 

* 

* 

* 

* 

« 

* 

* 

* 

« 

« 

« 

* 

* 

*  in 

*  a 

*  o 

*  x 

*  h- 

*  Ul 

*  X 

* 

•  a. 

•  uj 

•  *- 

*  in 

*  i 

*  -ft 

«  i 

*  x 

*  j 

*  z 

* 

*  I- 

*  *ft 

*  u 

*  “ft 

*  3 

*  a. 
«  x 

*  *-• 

# 

* 

*  « 


« 

« 

« 

* 

« 

* 

* 

« 

* 

* 

* 

# 

« 

* 

» 

3  # 

J  * 

ft  * 

* 

« 

* 

* 

* 

* 

# 

* 

* 

>  « 

I  * 

* 

:  * 

4  * 

« 

* 

* 

« 

« 

* 

# 

* 
ft 
ft 
ft 

ft  * 
ft  •"* 
ft  W  ^ 
ft  lO  +4 

ft  H  1 

ft  •  * 

ft  CNJ  — 

*  m  -» 

*  -*  * 
ft  •"* 
ft  »  — 
ft  CM -ft 
ft  in  0. 

5  -•  ii  to 

« 

ft  ~**~* 
#  -*  *  _ 

*  i  -«o 

ft  z 

*  —  cm 

*  x  x  o 

ft  -ft  4  10 


X 

4 

Ui 

ft 


4 

ft 

X 

ft 


3  — 
-  -J 
(0  ~ 

ft  <0 
ft 

3  -ft 
-  3 


H  4 

x  cm 
a  uj 
+  ♦ 


a 

UJ 
►- 
-ft  4 

•  -I 
-ft  3 
—  U 
H  J 
X  « 
0.  (-> 

•  uJ 
a 

O  4 


-ft  i/> 
—  ►— 
—  Jl 
X  3 
cl  m 
u> 
uj  ar 

k— 

3  J 
X  4 

X  z 
o*-» 
JU. 


a  ^ 

33 
4  * 
ft.  -ft 

*■* 

3  x 

*  4 
H  I 
»  ^ 

-ft  3 


4 

ft 

r 

♦ 


-ft  4) 


X  o 

♦  t- 

■J  o 

-  ID 


X 

^  UL  -ft 
•  II  ft- 
-ft  —  ►ft  - 

ini" 


Z  j£ 


X 
-ft  II 


II  II 
-ft  3 

*  S 

m  n 


•  » 


o  CM 

in  in 


o  o 
3  O 

CM 

in 


Z  CM  o  Z 

■“X  ^ 

ft«i  •  tJ 
II  II  <3  II 
*  —  ui  *-• 

*o 

n  o 

-4  —  *ft 

CM  — 

o  x  u.  o 

3  4  3 

3  in 

in  n 


-ft  CM  •  “Z. 

II  43  -ft  3 

3  “ft  —  * 


X 

4 

II 


fft 

lO 


o 

o 


rO  M  ^ 


II  II 

^  H 

-ft  -ft  x 
—  •>  -x 

<<  -ft  -ft  4- 

-  ■“  •'ft 

►_  lO  CM  -ft 


I  X  * 
z  x  Q-  to 


z 

CM  * 
*  X 


-ft  O  X 

3  3  X 


O 

o 


X 

X 


CM 

43 


X  ♦  ® 
V.  UJ  X  UJ 
N  3  UJ  • 
uJ  Z  I- 
►ft  in  *- 
--  *  f 
r  it 

O  CM  I 
J  * 

-ft  o  o 

45  3  fft 


-ft  10 

II  *■" 

*  z 

Z  X 

•  if  (0 

-ft  4) 

II  -«  -J 

►ft  —I 

O  4 

ioao 

o 


II  -ft  CM  X 
— 1  •  •  CL 
3  3  11 

m  uJ  uj  -ft 

vO  •  •  -ft 

-ft  4  Z  - 


o  X  X 
a  *-•  *-• 


lO 


V 

* 


I 

4 

UJ 

• 


4 

X 

X 

3 

4 
I 


* 

•O 


z 

ii 


x 

X 

» 


CM 

X 

4 

• 

X 

x  » 


>- 

* 


X 

4 

UI 

* 


4 

X 

X 

3 

4 

Z  2  » 


-ft  UJ  1 

£3* 

H  vO 


ii  ~  2! 

*  Z  fft 
V  -ft 
4)  II 
0  -ft  O 


ii  ii  x 
-ft  z  z 
>- 

.o  JO  II 
4)  4)  -* 


.  O 

*  a 


o  o 

a  3 

in  3  in 
4)  vO  4> 


3  Z  O 
O  v  3 

43 

4) 


0  3  2 

J  3  >■ 

r-  ® 

4) 


O  3  O  O  U 


—  o— ss-sssssaassaasaas^asssSSSSS 


) 


i 

i 


A-40 


TR  5341 


ft  ft  ft  ft 

ft  ft  ft 

ft 

ft 

ft 

ft 

•  ft 

ft 

ft 

ft 

ft 

ft  ft 

ft 

ft 

ft 

ft 

ft  ft 

ft 

ft 

ft 

* 

ft  ft 

ft 

ft 

ft 

ft 

ft  ft 

ft 

ft 

ft 

ft 

ft  -  ft 

ft 

ft 

ft 

ft 

ft  Z  ft 

ft 

ft 

ft 

* 

*  O  ft 

ft 

ft 

ft 

ft 

ft  *4  ft 

ft 

ft 

ft 

ft 

ft  H  ft 

ft 

ft 

ft 

ft 

ft  >i  ft 

ft 

ft 

ft 

ft 

ft  U)  ft 

ft 

ft 

ft 

ft 

ft  O  ft 

ft 

ft 

ft 

ft 

ft  X  ft 

ft 

ft 

ft 

ft 

ft  Z  ft 

ft 

ft 

ft 

* 

ft  O  ft 

ft 

ft 

ft 

ft 

ft  O  ft 

ft 

ft 

ft 

* 

ft  UJ  ft 

ft 

ft 

ft 

** 

ft 

*  a  * 

* 

ft 

ft 

X 

* 

•  » 

ft 

ft 

ft 

<t 

ft 

*  o  ~  ft 

ft 

ft 

ft 

CM 

ft 

*  *i  ^  * 

ft 

ft 

ft 

W4 

ft 

ft  D  +  ft 

ft 

ft 

ft 

X 

ft 

ft  O  Z  ft 

ft 

ft 

ft 

X 

ft 

ft  •-<  *■'  ft 

ft  Z 

ft 

ft 

Ui 

ft 

ft  X  >-  ft 

ft  o 

ft 

ft 

ft 

ft  Id  ft 

ft  t-4 

ft 

ft 

3 

ft 

*  X  _l  ft 

*  UJ  K 

ft 

ft 

z 

ft 

ft  <t  ft 

*  tr  < 

ft 

ft 

*4 

ft 

ft  £  Z  ft 

a 

ft  -d  CM 

ft 

ft 

< 

ft 

ft  X  M  ft 

3 

ft  X  >-• 

ft 

ft 

• 

ft  O  X  ft 

*  _J 

ft 

ft 

z 

ft 

ft  X  ft 

O 

ft  if)  < 

ft 

ft 

o 

ft 

*  rt  t-  « 

3 

ft 

3 

*»-*-• 

ft 

ft 

o 

* 

*  Id  U»  ft 

*  X  1- 

ft 

ft 

ft 

ft  X  D  ft 

(M 

3 

S' 

*  <  •-• 

ft 

ft 

X 

ft 

ft  3  ft 

X 

ft  K  Z 

ft 

ft 

-< 

ft 

#  «  a  * 

>- 

►4 

s 

ft  if)  1-4 

ft 

ft 

CM 

ft 

ft  h-  <  ft 

« 

ft 

ft 

ft 

UJ 

ft 

—ft 

ft  ft 

— * 

X 

H 

ft  X  1/1 

ft 

ft 

ft 

H 

*  X  Q  ft 

— * 

< 

X 

ft  UJ  UJ 

ft 

ft 

» 

ft 

«* 

ft  O  Z  ft 

CM 

3 

*1-0 

ft 

ft 

ft 

*H 

ft  <  ft 

* 

UJ 

ft 

•  1/10 

ft 

CM 

ft 

X 

ft 

ft  Z  ft 

M 

♦ 

ft  1 

ft 

-4 

ft 

< 

ft 

E 

ft  O  -  ft 

w 

w 

3 

ft  CM  Z 

ft 

c\l 

ft 

>4 

* 

-t 

ft  ►*  ~  ft 

X 

< 

M 

* 

*  1  o 

ft 

ft 

X 

ft 

=0 

H  IjJ 

ft  1 - ft 

3 

1 

» 

H 

ft  »4  >4 

ft 

O 

ft 

X 

ft 

o 

• 

♦ 

ft  (_><-«  ft 

t- 

S 

*1—1— 

ft 

t— 

ft 

d 

ft 

CM 

•H 

— * 

ft  Z  —  ft 

jJ 

3 

» 

CM 

ft  -J  o 

ft 

Z 

O 

ft 

ft 

* 

*H 

*3h* 

r 

►- 

S 

X 

ft  3  XI 

ft 

a 

O 

ft 

ft 

ft 

X> 

X 

ft  X  *-  ft 

M 

O  — 

< 

*  s;  ui 

ft 

H 

• 

o  o 

ft 

Z 

ft 

o 

<t 

H 

ft  <  ft 

a 

z 

w 

a  »-< 

ft 

• 

ft 

ii 

o 

3 

ft 

14 

* 

CM 

*— 

w> 

ft  <  ft 

3 

» 

>- 

•  r 

E 

ft  X  o 

ft 

3 

II 

•  — ■ 

ft 

«x 

ft  "M 

X 

ft  id  ft 

h  ^ 

"J 

<  T 

O  Z  3 

Z  Z  1 

ft  <  z 

ft 

.9 

— * 

O  —4 

ft 

c- 

ft 

X 

-t 

ft  l/IH  ft 

• 

< 

II  II 

—  *  *■*> 

ft  XJ  *-< 

ft 

«h 

-4  3 

II 

ft 

2 

ft 

30  UJ  UJ 

ft  >-4  <  ft 

3 

^  H 

-H 

♦  •-* 

^  H  ^ 

HHH 

ft  z  z 

ft 

-4  S 

o 

■—  • 

ft 

3 

ft 

o 

II 

*  3  * 

'il  I  h  II 

O 

II 

—  II 

>-4  II  M 

II  II  — 

*H7 

ft 

II 

CM  *-  i3 

•-4  K 

ft 

(J 

ft 

CM 

II 

—ft 

ft  <  _|  ft 

• 

♦ 

z  •-« 

3  3  ►*  w 

*  3  * 

-t  3  Z 

UJ 

ft  J  -I 

ft 

ft—* 

• 

*  3 

ft 

ft 

—ft 

H 

ft  4  ft 

X 

• 

w 

»  S 

► 

3 

ft  z  o 

ft 

0-4  0 

-4  • 

ft 

X 

ft 

#•% 

ft  X  >  ft 

r- 

<  »o  o 

•n  o  * 

^  S  - 

>n  x)  ii 

•  9 

ft  O  d 

ft  o 

.3 

X 

II 

ft 

i 

ft 

H 

H 

*  -4  Ui  ft 

< 

r«- 

II 

r*- 

il  r-  s-  r-  *  c-  ~ 

X  -4 

•  z  o 

ft 

* 

»  UI 

ft 

J 

ft 

1 

ft  ft 

►-« 

►- 

J  -c 

—ft. 

*4 

^  H 

^  4-4 

HHH 

3  »- 

ft 

ft 

CM 

M 

M 

ft 

ft 

9 

z 

ft  ft 

w 

II 

_i 

M 

M  M 

r-  Z 

ft 

ft 

w 

—t  — 

ft 

ft 

X 

i 

•  ft  X 

r 

4  O 

w 

o 

—  o 

e  o  r 

3  3Z  UJ  O 

ft 

ft 

3 

- 9 

X  x 

ft 

ft 

'Jl 

<  CM 

3 

3  3  3  0  3 

3  J>>  XU 

ft  ft  ft 

ft 

3 

V- 

X  -• 

ft 

ft 

ft 

j 

CM  <0 


«_>  L>  L>  L» 


K) 


S  to  o 

r-  o 

-4  -4  CM 


uuu  u 


-4  O  cS 

j  »  r* 

■mm  a 

M 

o  _>  u 


I 

i 

l 


1 


loioflrtioioioioio^  ^  i  ^  ^  £  <o  <0  -o  4>  £ 


A-41 


■"“*  “  *fl*wuj(»(S«ws*!BMs»MM»j»  ■*«  ttt  tfMmMnrmgDi 


TR  5341 


«  « 

* 

# 

• 

ft 

ft  ft 

* 

# 

ft 

* 

# 

ft 

« 

« 

ft 

« 

* 

ft 

* 

* 

ft 

o 

# 

* 

ft 

4k 

* 

* 

ft 

cv 

* 

« 

ft 

ft— 

# 

» 

ft 

►-I 

« 

* 

ft 

X 

» 

* 

ft 

CL 

* 

# 

ft 

* 

* 

ft 

10 

# 

* 

ft 

z 

* 

* 

ft 

H 

* 

* 

ft 

4 

* 

# 

ft 

K 

* 

# 

ft 

Z 

# 

# 

ft 

O 

# 

« 

ft 

o 

» 

# 

ft 

* 

* 

ft 

4“4 

« 

# 

ft 

3  CNJ 

* 

* 

ft 

•  10 

* 

* 

ft 

•ft  CM 

* 

* 

ft 

* 

» 

* 

ft 

CM  CL 

* 

* 

ft 

—  O 

* 

* 

ft 

*->  o 

* 

* 

ft 

X  _l 

* 

* 

ft 

CL  _ 

» 

* 

ft 

X 

* 

* 

A 

ft 

k  ►■» 

4-4 

~w 

* 

» 

3 

ft 

1 

— k 

Z 

k 

*  l/l 

* 

— 

ft 

1 

O  O 

•> 

iD 

#  a 

• 

5C 

ft 

«  UJ 

t0 

k 

*  o 

# 

4-4 

*  in  -♦  t- 

4k 

sT 

*  I 

# 

H 

ft 

•ft 

—  < 

X 

X 

#  f- 

—* 

* 

a 

* 

f\J  ►"•  _) 

4 

4 

o 

*  uj  eg  * 

* 

* 

X  3 

K) 

k 

40 

*  X 

i 

# 

k 

*10.0 

UJ 

X 

<Vl 

• 

» 

CM 

* 

ft 

o 

-1 

4 

*  a. 

« 

* 

3 

* 

m 

ft 

O  10  4 

X 

ro 

o 

*  UJ 

« 

* 

m 

* 

ft) 

ft 

JZU 

4 

UJ 

h- 

#  i- 

* 

<M 

— «* 

ft— 

CM 

ft 

•-4 

CM 

k 

*  i/i  i 

* 

H 

— # 

ft 

IL  4  UJ 

Ul 

X 

o 

*  i 

< 

* 

o  o 

• 

•-* 

X  o 

ft 

0  1-2 

4k 

4 

o 

*  C\J 

* 

J-  <\l 

H 

a 

ft  ft- 

ft 

X  4 

X 

CM 

*  1 

* 

k 

*— 

«» 

— % 

ft 

X 

O 

< 

UJ 

4"4 

#  X 

i0 

» 

O  -ft 

x 

z 

3  o 

ft 

0  3  — 

UJ 

k 

o 

*  _»  z 

* 

004 

k 

* 

o 

ft 

>■* 

H 

k 

<VJ 

#  2 

H 

* 

OJ  V 

X 

1-1 

ft 

•T 

~  ft 

I 

X 

• 

• 

< 

« 

— * 

M 

o 

'X 

w 

-k 

ft 

■>  T 

k 

k 

J- 

#1-1- 

* 

H 

k 

a 

k- 

~9 

•ft 

-ft 

(M  O 

ft 

4 

4 

o 

*  ►-*  2 

» 

•ft  o 

K 

k 

* 

* 

I 

ft 

►H  V 

- - 

• 

#  O  O  * 

•  o 

• 

X  *ft  *■* 

H  < 

• 

ft 

Jj 

<g  u 

J 

#  ►-*  3 

* 

H-  CM  -4  V0 

■J 

II 

II 

II 

II  <3 

ft 

-4- 

H  -1 

•ft  CJ 

a 

X 

•  -J 

*  o 

It  O 

>  t-t  ~3 

* 

—ft 

Ui 

ft 

< 

<\l  -t 

4 

UJ 

*  a.  eg 

» 

• 

—ft 

—  <\J 

3 

• 

ft 

z 

X  X 

X 

Q 

*  x  x 

» 

H 

HIO0 

K) 

k 

ft 

— 

r 

I-  JX  J  H  * 

J  II  J  -  t  *  ~ - -  _l  <VJw* 

O  <  eg  4  U.  *  *U.U.X04000XU.* 

4)^IO-*  ##  «  3  4  -•  # 


*  # 


* 

* 

« 

* 

» 

» 

« 

« 

* 

* 

* 

* 

* 

* 

* 

« 

* 

* 

« 

4k 

M 

J- 

t-4 

ft 

z 

-}  — 

3 

« 

X 

-  3 

k 

ft 

+ 

jC 

H 

-  ^ 

k— 

« 

~i 

X  — 

X 

» 

4k 

4  X 

^  4 

* 

1-4 

rsj  4 

3  ft 

* 

UJ  jj 

-  O 

• 

X 

ft  ft 

►ft  Q 

* 

4 

“  O 

* 

U 

I  • 

» 

ft 

4  CM 

« 

0“+> 

t-4  H 

1  ft 

* 

C\J 

W 

—ft  —ft 

* 

' — 

X  X 

3  3 

* 

<1 

4  4 

k  k 

* 

X 

ft  ft 

H  •-« 

CL 

4-4.  4-4 

ft-  ft— 

* 

_J 

-ft  CM 

H  H 

* 

< 

X  X 

* 

♦ 

4  4 

i  -f 

• 

4-44 

X  X 

4— ft  —ft 

* 

~i 

CL  X 

3  3 

* 

• 

-J  J 

k  k 

• 

H 

4  4 

—I  »-• 

* 

44 

ft  ft 

k  k 

* 

X 

—ft  —ft 

*-♦  Al 

« 

4 

->  -J 

ft—  ft— 

• 

A» 

»  •> 

►-« 

* 

UJ 

M  »-H 

X  X 

♦ 

ft 

k>  k 

4_! 

X  CL 

* 

~4 

^  (M 

• 

II  II 

« 

w 

H 

— k  — k 

* 

4-4 

M  H 

II 

O 

* 

< 

r  r 

1 

k  k 

* 

X  2  X  X  Z  .X 

X  CL 

Z 

H  H 

* 

k  k  J.  4k  k  4k 

II  II 

-  O' 

k  k 

* 

H  H  J  H  H  H 

—ft  —ft 

H  H 

CM 

II  II  4  II  II  II 

II  ^1 

k—  ft— 

» 

-•  X  II  -•  aC 

k  k 

H»  H 

* 

— k 

M  H 

O  X  X 

* 

fO  ro  0  r>  40  r> 

k  k 

S3  3 

X  X 

* 

-ft  M  »  H  H  H 

4-ft  M  4ft 

ft  tM  r\J  -4  fg  c\j  rg 

ft- 

eg 

ft  — 

ft  O  D  -«  O  O  O 

M  H 

x  r  o 

ft 

3  X  X  o  X  X)  CL  a 

GO 

eg 

•ft  —4 

4-1 

n 

o 

»ft 

10  O 

o 

o 

<A| 

eg 

u  u  u  u 

eg  g 
g 

eg 

eg 

-o 

K' 

in 

o* 

o 

-* 

-♦ 

-4 

ig 

g 

eg 

g 

U  3  U  3  <3 


H  -I  ^4  H  H  Sl^M 


A-42 


i03  218  PI ( 1 # 1 ) =0 .00 

*,54  uu  £20  K:l»KSUP 


TR  5341 


ft 

»  » 

» 

• 

ft 

ft 

3 

ft 

ft 

. 

ft 

ft 

-H 

ft 

ft 

— . 

ft 

ft 

X 

ft 

ft 

ft 

ft 

<4 

ft 

ft 

UJ 

ft 

ft 

ft 

ft 

ft 

^*«» 

ft 

ft 

cH 

ft 

ft 

w 

ft 

ft 

4 

ft 

ft 

X 

ft 

ft 

X 

ft 

ft 

3 

ft 

ft 

< 

ft 

ft 

1 

ft 

ft 

ft 

ft 

3 

ft 

ft 

4fc 

ft 

ft 

(4 

ft  1/1 

-~v 

ft 

♦  Or) 

ft 

>- 

ft  o 

1 

ft 

ft 

ft  I 

ft 

ft  h* 

ft 

ft 

3 

ft  II) 

ft 

ft 

. 

ft  X 

ft 

*-i 

ft 

X 

ft 

ft  X 

< 

ft 

X 

ft  XI 

ft 

< 

ft  ►— 

ft 

3 

ft  v/1 

n 

ft 

**) 

ft 

ft  1 

2 

ft 

•-4 

*■“ ft 

*-* 

ft  14 

ft 

E 

3 

4 

ft  1 

4 

ft 

H 

ft 

w 

ft  *-• 

1- 

ft 

w 

t 

< 

ft  h- 

~r 

ft 

I 

3 

>- 

o 

X 

*  J 

6 

ft 

4 

# 

W» 

ft 

O 

X 

*30 

ft 

ft 

◦ 

— » 

• 

J 

ft  X 

ft 

ft 

3 

o 

4 

-ft 

ft 

m 

ft 

3- 

iH  iH 

< 

. 

II 

1 

«• 

ft  X 

E 

ft 

00 

H 

<fc 

3 

H 

»■» 

2 

ft  4 

4 

« 

ru 

-H  H 

w 

t— 

H 

►ft 

•  tt 

ft  UJ 

ft 

.4 

I 

* 

»H 

3i 

*  D 

3 

ft  2 

Cft 

ft 

◦  vO 

2  C4 

<1 

* 

* 

X 

♦  HS 

ft 

y-  <4  <i  i 

-f 

< 

►ft  3 

ft  3 

4) 

ft 

•ft 

♦  < 

>- 

1 

•H 

S 

I  X 

% 

ft  2 

<4 

ft 

◦ 

AJ 

^  w 

1 

* 

»H 

a. 

►» 

X 

C4 

ft  ◦ 

ft 

◦  4> 

H  \ 

* 

. 

r 

H 

►4 

ft  — 

ft«^ 

ft  2 

i 

ft 

<4 

^  O 

M 

a 

. 

— * 

X 

—  □ 

H-4 

ft 

O 

ft 

31 

^  a 

T*  ~~9 
ft.  •_ 

2  ♦ 

2  2  — i  2 

o 

2  X 

-ft  it 

2  I 

ft  1- 

3 

ft 

w  O 

j  :  z  ^ 

il 


<s>  -« 

—  il 

2  H 

u. 

o 

31 
J  31 

3 

-t  O 
O  3 


II  •«' 
■J  2 
>- 
-t  II 
CM  - 
N  -• 


II  II 
Id  -•  3 
◦ 

2  0  0 
►*  <o  m 

t-  4  ,3J 


H  -H  *  H 

a  II  r?  II  -l 
II  -•  'd  3  » 

•  H 
-•  4  fO  - 

-•o  —  os- 

-♦  rsl  -t  «\|  — 


»  * 

l  3  ►« 
Ui  II 


a 
—  n 

3  — 


►  a.  f  WJ  # 

^  II  #  O  X  ft 

II  —  O  ft  ►-<  # 

TTMIJl-t 
-  iH  ft  X  X  * 


H-  O 
<  r* 

h  3J 


o  ◦  ◦  o  -i 
3  -  O  3  3  X 


3U.O 
3  -•  3 


«-»  O 

CM  v\| 

4  4 


o  in 

■o  -o 

4  4 


I  —  2  a.  ◦ 
a.  >■  >■  h  j 

->» 

K)  (V 

K)  »») 

4  4 


O  ►» 

7>  —  O 
Mi¬ 
ll 

Oi\IO 
3  3  <3 

O 

3 

4 


ft  5l  Ld  * 

ft  2  * 
*  * 
*  * 
ft  ft 

ft  ft  ft  ft 


CMI  •  »  *  »X  » 

•  n  ^  H  *H  H  H  <1  *"4 

I-  (M  II  II  rO  II  II  II  II  II 

xi  ◦ 

3  •  —  3 

2  h  *  •  ;r  j-  3-  »  r- 

-ft  J)  I  ftft  )i(M)  flft'O 

>-  -«  2  r\|"\J  vl  —  t\i 

2  —  —  CM  <)  iM 

◦  xu.  e  r  a  o  ◦  ◦  r  o 

3  J  4  3 

OJ  3 

O  3 

V  4 


◦  -ft  — i  X  <  2>  3 


O 

0 

4 


U  J  uu 


CM 

vC 

4 

4 


X)N4)jro»ft.M'OJ-inft)r-  o^o-ft^ioj-noi^-o 
OOOOO-ft-ft-ft-ft-ft-ft-ft-C-ft-trMCSIMaiM  4  31  3)  (4  .  4  *>  «3  iO  •)  lO  fl  O  lO  4 
4  (4  4  ll  >\l  M  M  M  Al  4  4  -4  31  M  4  V  4  4  4  4  4  4  4  C4  3J  4  4  4  4  4  4  4  4  4 


A-43 


TH  5341 


#  #  * 
* 

» 

» 

* 

« 

* 

» 

* 

ft  <M 
ft  *• 
ft  -H 
ft  * 
ft  O 
ft  II 
ft  _l 


« 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 


ft  *  4)  ft 


-5  ~  3  ~  3  ~ 
-  3  -  3  *3 


X-I—X*- 
<  X  «  X  <  X 
f\l  <  IM  <  <M  < 
txJ  Itl  UJ  U  aJ  LU 

ft  ft  ft  ft  ft  ft 


ft  -ft  ® 

* 

«"k  — X  ^k 

*  _i  g 

* 

333333 

*  » 

« 

k 

k 

k  k  k  k 

ft  g  cr 

* 

H 

H 

I 

I 

I 

I 

ft  ft-O 

ft 

k 

k 

k  k  k  k 

*  -4 

* 

.4  g  -«  g  -*  cm 

*  x  g 

* 

w 

w  w  w  w 

*  a  co 

* 

^k 

rtHIMNPOlO 

*  g 

« 

3 

3 

* 

x  x  x  x 

*  i/> 

» 

k. 

ft 

* 

«  *  «  • 

#  z  0. 

* 

H 

•4 

*k 

#«k  *k  ^  *^k 

ft  M  0 

* 

w 

k* 

rlftlHNHN 

#  <  0 

* 

g 

g 

W 

w 

•w  w  w*  «-* 

ft  K  _l 

« 

r 

X 

< 

4 

4  4  4  4 

-k 

-k 

ft  Z 

* 

<  — 

< 

X 

X 

X  I  X  X 

3 

3  3 

ft  O  X 

« 

4-  — 

+ 

cl  a  d  a.  a.  a 

k 

k  k 

3 

ft  O  •-* 

« 

3  3 

333333 

H 

M  — « 

k 

ft 

* 

333- 

k  k 

3 

<<<<<< 

k 

k  k 

* 

#  ~  a 

* 

»  »  *  H 

H  H 

k 

+ 

ft 

ft  ft  ft  ft 

ro  r>)  K) 

ft—* 

ft  3  UJ 

* 

H  H  H  w 

w  w 

M 

X 

*  -i- 

« 

—  -  —  X 

g  x 

w 

3 

3  3 

H 

g  to 

< 

*  p-«  < 

* 

X  X  X  « 

1  < 

r 

k 

k  k 

X 

X  X 

ft 

*  -  _i 

« 

<441 

«i  ♦ 

4 

►4 

►4  ►H 

ft 

4-  + 

*JO 

« 

ft  ft  ft  — 

■ft  — 

ft 

k 

k  k 

— * 

*k  -k. 

* 

#  ~o 

• 

0003 

—  3 

0 

H 

g  <o 

3  3  3 

k 

ft  -«  _i 

» 

3  3  3  * 

3  * 

0 

w» 

k 

k  k 

M 

^  .n 

#  1  < 

« 

0  in  iflM 

•  H 

in 

H 

►4  H 

H 

1-4  H 

ftp*  <*ft 

r-  x>  *  a.  0 

« 

•  •  w 

H  w 

• 

X 

X  X 

k 

k  k 

g  to 

g  g 

# 

* 

-4  0  0  1- 

w  t— 

-4 

X  3 

LS  3 

^4 

g  0 

I  X 

ft  *  UJ 

* 

1  1  ft  —1 

ft-  —p 

ft 

•  II 

II  II 

w 

«x  < 

0  0 

*  in  or 

« 

**k  «*k  -k  ^ 

-•  z 

#*•* 

*-■» 

-4  — 

^k 

H 

H  1-4 

ft  ft 

»-  ft- 

*  < 

* 

3  3  3  3 

Z  3  3  3  3 

II  3 

3  3 

1 

X  X 

-ft  2 

*  g 

» 

k  k  k  W 

J  W 

k 

* 

k 

Z  3  * 

k  k 

3  3 

.3  3 

■)  ft 

00#  — 

-•  -ft  -«  * 

ft  ft 

►4 

*4 

►4 

*  ft4 

►4  H 

3 

11 

11  11 

»  ^ 

0  0 

#0.-1 

# 

—  —  —  O 

O  O 

k 

k 

- 

k  k 

Z 

4*k 

*k 

—  I 

*  0  ♦ 

# 

ft-  ft-  H  O 

3  0  g  10  ft- 

II  ft-  -4 

g  po 

H 

3 

33 

— *  ^ 

— k 

*  0  z 

* 

H  H  H  • 

•  #• 

k^ 

►4 

3  g  — 

w  w 

>- 

k 

k  k 

4  ft* 

-4  O 

*  _i  — 

« 

zz. 

4 

z  z  z  g 

NAIHH 

z 

►4 

N4  1-4 

z 

M 

*-4  H 

x  ft- 

*  >- 

« 

* 

* 

3  3  3  1 

1  t 

X 

X 

5 

k 

•0  0  X 

x  1 

0 

k 

k  k 

<  cr 

•  • 

*  u. 

* 

-« 

II  II  II  II 

II  II 

11 

11 

11 

-4  ft-  3  O 

3  3 

3  -# 

g  10 

II  UJ 

1-  a 

ft  0  _l 

• 

II 

;i 

— *  <*k  — «» 

«*k 

— % 

11  g 

— * 

•W  k# 

-ft  :»  UJ  ID  Id  ft  < 

ft  «-*  3  3  3  3  3 

3  3  3  3  3-* 

H 

►4  *4 

3  Z 

3 

•  • 

*  a  z 

ft 

k  k  k  k 

*  * 

k 

k 

k 

0 

X 

X  r 

•  *"4 

T 

*  z  -* 

ft  0 

T 

T 

T 

I 

0 

"4  H 

>-4 

*-| 

-» g  3 

a 

>  3 

1-4 

—•in  ft- 

ft  UJ  u. 

*  i-  ft- 

k  k  k  k 

k  * 

k 

•ft 

k 

ft- 

—  _J 

»- 

>-•  »-4 

ft 

*  g 

g-rgo-ftgo-* 

g  <o  g 

4  J  Z 

ft—*  k* 

« 

ft 

w  w  w 

4-4 

w 

w 

X  < 

ou.ii. 

♦ 

* 

0 

0  •*  —  g 

g  g  po  0 

ro 

0 

<  O 

u  —•  »-* 

»  «  « 

* 

3 

3 

<  #  c  # 

s  « 

* 

X 

3 

r* 

(M 


tO  4- 
X)  JO 

M  g 


o 

ft- 

g 


4 

ft- 

Ol 


u  u  u  u 


3»o^4.ag«04mg)N«9«o*ftcgK)4i><0r>(0a«O'ftcgi04inv0r^a09>O'ftAi04 
'a4444444444in.ninin.n;n(n;oinin.og)<0(0«NO<og>g>«or‘'f>»f‘«i>*r* 
g  \i  g  g  g  g  g  g  vg  g  g  g  g  g  g  g  g  g  g  g  g  g  g  g  g  g  g  g  \i  g  g  g  g  g  g  g 


A-44 


t7b  «l73  CONTINUL 

*7to  c’U  Cwi,T  1NUE. 

<.11  UO  ^7t>  L=i#i 


i 


TR  5341 


«  *  4  * 

M 

ft  ft 

•  ft 

X 

ft  ft 

4 

ft  ft 

lil 

7 

ft  ft 

* 

% 

ft  ft 

<4 

ft  Z.  ft 

C4 

ft  0  ft 

>- 

ft  Ht  ft 

4 

* 

ft  1-  ft 

ft  *-*  ft 

a 

7 

ft  l/>  ft 

_l 

ft  O  ft 

X 

H 

ft  X  ft 

1 

ft  X  ft 

X 

ft  0  ft 

7 

4 

ft  0  ft 

Id 

ft  Id  ft 

* 

ft  0  ft 

w 

A 

ft  ft 

>- 

(4 

ft  u  ft 

ft 

•w> 

ft  *-«  -t  ft 

4 

ft  0  +  ft 

■7 

I 

ft  0  X  ft 

«• 

X 

ft  »-4  —  ft 

-> 

t— « 

-1 

*  ac  >-  # 

v 

4 

ft  UI  « 

X 

1 

ft  X  -1  ft 

4 

ft  4  ft 

71 

7 

ft  X  Z  ft 

w* 

Ui 

» 

•  ac  h«  « 

•-ft 

* 

•-< 

*  0  lL  ft 

4- 

ft  lL  « 

3 

-4 

>» 

*  a:  >-  * 

# 

* 

*  u  1/1  » 

4 

** 

ft  0.  7  * 

X  x 

7 

*  7  * 

•> 

*  X 

* 

*  ►  a  * 

M 

_» 

H 

ft  H  4  ft 

w 

7  4 

w 

ft  ft 

*o 

—  I 

X 

ft  U  Q  « 

X 

7  — 

4 

ft  O  Z  * 

< 

*  >-< 

(4 

ft  4  ft 

1 

< 

Id 

ft  Z  * 

7  X 

ft  <2 

ft  O  •  * 

-  tr 


-  H  II 

-4  7 

ft  ft-  ~ 

ft 

H 

*  •  — • 

—  *- 

ft  U  M 

ft 

»  jC  —* 

4  Ul 

ft  Z  ~- 

» 

j 

>-  M  w 

x  t: 

ft  7  K 

ft 

-  M  Z 

X 

ft  U 

ft 

M 

Z  X  >- 

J  -> 

ft  4 

ft 

X 

•  Q. 

4  O 

ft  4 

ft 

X 

:nx  ’  ♦  — 

Z  1 

ft  UJ 

ft 

II 

*»>•««  lO 

ft  M  • 

ft  U)  1- 

ft 

H  4(5  H  M 

*4  ►"«  v3l 

*►4  4 

ft 

4 


7 

» 

tO 

>- 

ft 

~>  ~ 

*  7 

<4  » 

•—  t-t 
>-  — 

*  4 

~  I 

7  _j 

--4 

4  K 
I  4 

3  ♦ 


II  II  ••  ”>  II  II  «—  II  «—  • 

H  T  ^  «Hi2T23 

►4  U.  >-  Hi 

^in/)»DOOiMii  • 
JHO  CO 

4  i4  4  —  i4  4  J  C4  <-•  Z 


(7 


w  W  # 

£ 

— 

7 

7I7740Zu.~ 

O  7 

7 

7  Z  a. 

7 

7  X  7  7  'J  7 

►  Hi 

7  3 

7 

□  >-►-« 

* 

lO  J> 

C4 

O 

m 

4) 

ft*  r> 

O 

OJ 

D 

CO  > 

(4  4 

M 

4 

4 

c4  4 

—  Klllll'-'UJ*  7  *  x  *Z*1 
")  1  J  O'  -4  7  2  •  ft  4  -I  ft  ft  2  2 
iO<4  >-  *  «  *  ~XU.lL 

."Mi:  UJvOIIP-ftlL^ft-l+44 

^HO(OS''<  t  h  ii  #  «4H 

ft  t  H  M  <  4 

*  *  ft  ft  ►-  I-  O  J 

O' 

M 

uu  uu 


’  o  '  0.  -4  S  O 

•  D  •  O 

H  .  H  +  H  H  • 

II  OC>  II  ^O.  II  O 

•-«  Q  II  O  ^  II  H  II 

?MI  »  O'  II  »  0*  » 

v4  ,~i  (V - 4  4 

3  0  0.  7  3  a.  7  CL 

4- 

cr 

4 


9  9>OH 
P-  O  JO 
4  J4  .4  4 

* 

\ 


M’O^lfl^NOO'O'^N'O^ 

JT  4)  D  45 

4  4  4  -4  4  4  4  M  4  *4  4  4  4  4 


flPBJOrt.ilflJflJMOJ'O 
3»03tff»0000000000-4 
\l  4  4'4  0<lf0  0^>00,0-0'0  7 


/ 


A-45 


( 


TR  5341 


Q. 

# 


IO 


X 

0. 

♦ 


CO 

ft 

—  CO 
3  O' 

-  CM 

H 

-  o 

CM  t- 

n  O 
X  O 
Ql 

+  «** 
to 


z  CM  G» 

*  —  UJ 
• 

ii  a 
->  ii  a. 

—  ui 
m  *+  k 
o*  *  in 
CM  CM  X 

OHb. 
Ql  •* 


CM 


0. 

ft 

X 


CM 


<L 

ft 

X 

Z  I 


II  *•* 

»•  X 
>- 

0  II  Z 

CJV  X 
(\l  h  Q 

O  2  ul 
3  »-  uC 


II 

M  z 
>- 
O'  II 

cr  — 
CM 

o  z 

L»  >* 


*  ft 

ft  ft  ft  * 

1 

1 

ft 

ft 

ft 

ft 

• 

ft 

ft 

I 

\ 

ft 

ft 

J 

ft 

rO  ft 

ft 

1  ft 

I 

ft 

—  ft 

! 

* 

ft  ft 

ft 

ft  ft 

I 

* 

—  ft 

4 

ft 

X  * 

j 

ft 

<t  ft 

ft 

—  ft 

\ 

ft 

1  # 

i 

ft 

•  1  * 

i 

i 

ft 

•  to  ft 

ft 

•  X  * 

\ 

ft 

o  <  * 

\ 

ft 

Ul  »  * 

3 

ft 

a  —  * 

1 

ft 

OX* 

ft 

i-  <  * 

i 

ft 

Ul  IO  ft 

j 

ft 

«•  ft 

i 

ft 

Z  UJ  CL  * 

* 

ft 

02  X  * 

ft 

h4UJ* 

ft  UJ  I- 

*  X  <  in 

ft  uJ  CM  *- 
*XH  J 
ft 


I  # 
1  * 
X  * 
<  * 
IO  * 


ft  i/I  <  Ul  Ui  * 

*  l-  M  J  »  * 

t  1  h  OC  ^  * 

*  I-  Z.  &  <£  * 

ft  Ul  >-«  Z  CM  * 
ft 


0.  in  M  X 
Jlil  O  X 
ft-  O  -I  Ui 

IflQ  J  I 
I  O  I 
K1ZU.I 
I  o  < 
Jj  CM 


Ui  * 
3  * 
Z  Z  * 
X  *-•  * 
31-  * 
k-  Z  * 
ui  O  ft 
X  <0  * 


ft 
ft 
ft 
ft 
ft 
ft 
ft 
ft 
ft 

ft  _J  o 
ft  3  Ui 
ft  X  l/) 

ft 

ft  AC  i9 

ft  <  z 

ft  Ui  *■« 
ft  2  Z 
ft 


ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 


ht-IUJ* 


r-t 

~  -4  X 

+ 

w-i  »  o  C 

CL 

Ul 

*  -t  in  * 

•u  •~  ,0  ■" 

ft- 

U) 

H 

CM 

„  -  IH  —* 

X  X  <  *  ?  ~ 

II 

a. 

3 

X 

■O 

*  *t  lii  ft* 

CM 
O  o 
ft-  to 

•»  UJ  ft  — 

H  *  -IO(M 

«.  —  *4  «t  O  x 

•4  -«  -  *  < 

3 

*: 

X 


X 

ft 


J  H 
Z  19 
O  Ui 
Z  B 


-  X 
to  «y 

O  — 
fO  Ol 
X 
I  ii  ui 

X  I 

o  » 

Ul  X 
±1  < 
O  ‘U 


ft  ft  ft  ft 


a  z 

Ui  * 

•  •M  O 


3  CM  — 
ID  O  < 
IO  — 

a. 

H  »  x 

CM  Ui 

O  3 

(O  II 


—  I 

I  « 
<t  CM 
Ui  Ui 

II  II 


X  HO 
UI  X  »-«  •  o  *"* 
ft-  lii  O 


.  Ul  o 
l/IQO  II  r 
*  ZCM  '-Ul  I 

II  >ftlO  hhZ  - 

X  —  —  —  O'  X 

3  U.  O  Z  U.  U.  < 

J>-HHil 


—  f.9 

X  Ui 
<  • 
II 

*■»  X 
•ft  Ui 

o 

•ft  z 


o 

a 

o 


II  to 
-  o 

•ft  K> 


x  r 

4  < 

cm  to 

lii  ui 


CM  — ' 
X  u. 

<  >-ft 


to 

X  o 

<  o 


Z  Z  Z  cm  .2 

»  X  - 

H  H  rft  4  ft 
II  II  II  II  II 
►ft  3  X 

3 

O  O  O  *  tO 
■o  IO  O  ft  o 
(OO  iO  -to 
‘M 

O  O  O  I  o 
3  3  0  4  0 


o 

. :  Q 

■h  • 

•ft  o 
ii  ti 
3  — 
3 
O  - 
to  1-ft 
tO  — 

s 
3  X 
3  4 


40 

M> 

<S  O' 

o 

o 

CM 

CM 

o 

<7* 

O' 

O'  O' 

o 

tO 

K> 

CM 

CM  M 

fO 

u  L>  O  u  O  U 

O 

CM 

o 

iO 


o 

to 

ro 


sSKSSSSamSSSHSHHSMSIHSHH* 


A-46 


TR  5341 


* 

* 

* 

* 

* 

# 

* 

* 

ft 

« 

* 

« 

* 

« 

ft 

* 

« 

4a  4a 

4a  4a 

^a  4a 

* 

* 

33 

-33 

4a 

3  3 

4a 

• 

ft 

a  a 

3  -  - 

3 

a  a 

3 

* 

* 

*  x 

a 

X  X 

a 

# 

* 

<4 

■*: 

V-4 

X 

* 

-a 

* 

I  X 

—  XI 

X  I 

v> 

* 

3 

ft 

4  4 

X  4  4  X 

4  4  I 

* 

a 

* 

-a 

K»  CM 

t  K>  lM 

4 

3  CM 

4 

# 

K-i 

* 

3 

UJ  UJ 

UJ  UJ  U  UJ 

UJ  UJ  UJ 

* 

a 

« 

a 

ft  ft 

ft  ft  ft 

* 

ft  ft 

ft 

* 

* 

M 

4a  4a 

^  4a  4a 

4a 

4a  4a 

4a 

* 

* 

*  x 

z 

X  X 

X 

* 

M 

* 

OJ 

—a 

a  a 

a  a  a 

a 

a  a 

a 

# 

X 

* 

r 

3 

H  M 

MHH 

M 

M  H 

M 

* 

a. 

ft 

4 

a  a 

a  a  a 

a 

a  a 

a 

« 

ft 

ft 

H 

*-t  CM 

IOHWK) 

1 

2 

3 

* 

X 

* 

o 

W 

W  »4 

w  V*  V-4 

a4 

ft-4  V4 

S4 

* 

H 

* 

^  Cj 

CM 

•H  «H 

<rft  Cft  CM 

CM 

IO  K)  IO 

• 

i/> 

* 

3  3 

T 

*  rs 

Jt 

*  ^ 

■si 

* 

* 

»  »  K) 

4 

*  * 

ft  ft  ft 

ft 

ft  ft 

ft 

# 

t- 

» 

H  H  | 

ft 

4a  *a 

4a  4a  4a 

a 

4a  4a 

4a 

4a 

* 

_i 

* 

w  w 

o 

CM 

1*1  •-*  CM 

K> 

•4  CM 

m 

z 

* 

3 

* 

CM  CM  3 

O 

w  v» 

W  -4  W 

V-4 

w 

a 

* 

l/l 

« 

r  i 

• 

4  4 

4  4  4 

4 

4  4 

4 

u> 

* 

UJ 

« 

4  4  —  —  i-« 

to 

X  X 

XXI 

I 

X  X 

X 

a 

* 

cr 

ft 

4  ft-  —  ~ 

4 

Q.  CL 

q  a  c. 

a 

x  a 

X 

3 

4a 

*a 

* 

ft 

- >  3  3  3  X - 

4-* 

_J  _l 

_i  _j  _ i 

_i 

-J  _i 

_l 

a 

«-ft 

ft— 

H 

* 

a 

* 

33  13  -  -  -  433 

3 

4  4 

4  4  4 

4 

4  4 

X 

X 

a 

a 

* 

<\j 

ft 

a  a  a  a  H  M  M  t  a  a 

a 

4  4 

4  4  4 

4 

4  4 

4 

+* 

a 

o 

H 

* 

a 

« 

»-• 

4a 

4ft 

4a 

X 

O 

a 

a 

* 

H 

* 

—  —  —  —  X  CM  X  O  ’ 

V4 

1 

3 

3 

4 

a 

a- 

* 

a 

» 

IXXI4X4  -XX 

X 

a 

a 

a 

ft 

X 

i 

W 

* 

o 

* 

4444  1  44*44 

4 

H 

*-l 

M 

4a 

4 

4 

* 

II 

* 

ft*  ft  ft  —  4  ~  1  ft* 

ft 

a 

a 

a 

X  ID 

w 

at 

* 

» 

x 

* 

00003  —  3—  00 

O 

•H 

CM 

m 

«* 

iJ 

i 

a 

# 

* 

0300  -  1  -130 

Q 

w 

w 

H 

a 

4 

;t  o  * 

a 

« 

J1  J)  ill  JlH  -  ft  •  t f)  f> 

in 

W 

»-• 

M 

o 

—  X 

K> 

X  f* 

# 

—a 

* 

• 

l 

X 

l 

Al 

*— 

fO  4 

UJ 

4 

O 

* 

x 

* 

'-COO-tK—  t—  —  oo 

CM 

•ft 

3 

J 

3 

X 

I  CM 

a 

a 

* 

a 

» 

1  1  4  4  ft  1-  »-•  t-  1  4 

4 

a 

II 

II 

II 

4 

O  4  UJ 

I 

r  o 

*  to 

ft 

—  —  —  —  Zt-tZ*-*  —  —  — 

4a 

H 

4ft 

4a 

4  <!> 

a  a 

4 

4  4- 

* 

w 

ft 

3  3  3  3  3  Z  3X3333 

ii 

3 

3 

3 

4a 

X 

i\J 

IO 

« 

M 

ft 

»  -  —  —  —  3  —  3  -  —  — 

a 

T 

x 

a 

a 

a 

3 

«-a 

••  < 

UJ 

UJ  o 

» 

X 

ft 

HH  HM  *  ft  ft  ft  l-t  -t  ►-* 

>4 

a 

H 

H 

a 

*  UJ 

a 

a 

© 

« 

0L 

ft 

wwwQOOO-''  - 

H  O' 

a 

a 

a 

M 

X  - 

OJ 

tO 

« 

ft 

>- t—  b-b-oonah-b-^h- 

ii 

o  -c 

CM 

V4- 

• 

4  X 

X 

X 

4a 

#  UJ 

ft 

H  H  M  •  •  •  »H  H  ^ 

H 

3  n 

t4 

ft4 

ft-4 

x 

3 

W  a 

a 

a 

o 

*  f- 

ft  Z 

z 

Z  Z  .Z  Z  CM  CM  -M  CM  :r  T  *4 

7 

•9 

»-« 

w 

M 

* 

IUU-< 

4 

4 

* 

4 

ft 

a 

a 

3  3  3  3  1  1  1  13  3* 

3 

a 

<o  O  X 

I 

X 

4 

• 

ir  — 

—4 

V4 

a 

• 

-1 

ft 

II  II  II  II  II  II  II  II  II  II  II 

II 

-t  o 

O  3 

3 

3 

II 

n 

UJ  uJ 

UJ 

UJ  C5 

« 

3 

ft 

ii 

ii 

il  m 

^a 

X 

>  Q 

a 

O  UJ  ft 

O 

ft  i-i 

•5333333333333ft 

i  J 

Z  4 

4 

X  4 

•  « 

3 

ft 

a 

o 

!*)  *•_)*-«  i_L  I  D.  ♦  Q.  t  < 

-O'-*  Z  4  CM  ►-  ft  u  ft  o  o  ------  ------  o 


JI  JH# 

^  -  JJII  J  llj'-f  #  - ■ - - - - - - - - 

3X*-*OL>IOX<J>-'***3:j.«*.««iSft.*:5*ft**a 

#  *  ft  ft  *  ft 


IO 

to 

&  o 

to 

o 

O  vO 

>o 

o 

o  o 

o  o  o 


r^t><y*o-^vo^-in<oh-®^o^.'\jio*/>yOr^x)(7«o^cMK)4-in<or*-.oa‘o^cM'0^tn 
j-j-4-^)in^j)inin/)inx)X)43^s04)  <o.o>o>o<OM>r-rvr-p-f-r-rvfvfvrv:o.o.o  a  o  o 
OOOO-OO'OOlO'OiO'OO'OOOOO'OOtO’OiO'OOrOOiO-OO/OO  o  o  o  -o  o  V) 


A-47 


TR  5341 


3  3  ~ 
-  •  3 

ac  x  * 

^  ^ 

x  i  »— 

5  **  x 

bJ  bJ  UJ 

*  *  ft 

*  *  * 


< 
X 
X 
O'  -I 
o< 
w  ft 

O  3 


o  «. 
13  * 


I 
C3 
•  II 
UJ  — 
2  ■) 


a  • 

*  « 
H 

U.  I 

*+  3 


l  Ol  fO 

i  * 

•  * 

tv  <n 

<  < 

X  X 
0.  0. 
-J  -I 
<  < 
ft  ft 


;* 

Tt 

ft 


3  3  3  -. 

•  •*  4- 


3  3 


*  it 


X  X 
a  o 

*  # 


»o  * 
I  X 
<  «* 
I  I 

T  -3 


* 

ft 

ft 

ft 

• 

* 

* 

* 

* 

« 

* 

* 

* 

* 

* 

* 

ft 


» 

* 

* 


*  » 
* 
ft 
# 
* 
ft 
* 


X 

< 

<v 

UJ 

* 


« 

* 

* 

* 

» 

« 

* 

* 

* 


«  ft 


*  I 

•  - *-  3 

r?  o»  w  ii 

♦  ♦  ♦  3 
333h 

^  ^  Oh  % 

H  W  M  ^ 

>  •  K  V> 

H  V  fO  H 

- - X 

MHHO 
XXI 
U  3  3  'J  ^ 

^  »l  II  II  * 

4~  *m*‘  • 

«"3->->o  a 

?r  *  -  -  ui  uj  3 

Z  hm  w  v  ?  -■ 

P  *  *  •  0.  2  *  “t  t 

m,**  ,  y  (*  ii 

*  X  X  XOO  * 

33  3hu^>  00o 
*-•  -*  *4  -» 
►“  »0  *)  X) 

-•W 

o  O  O  o 
o  a  qq 


* 

tv 

* 

# 

< 

# 

I 

« 

X 

* 

-J 

* 

< 

* 

1 

* 

# 

3 

# 

* 

* 

K> 

ft 

— 

ft 

>- 

ft 

ft 

ft 

ft 

3 

ft 

ft 

H 

ft 

N-# 

ft 

X 

ft 

«T 

ft 

bJ  — 

ft 

*  3 

ft 

^  % 

ft 

ro  ^ 

ft 

w> 

jt  ^ 


X  X 

x  a 
ii  ii 
— % 
"3  -> 


5£  * 


X  X 

x  a 


O  r-i 
2  •  . 
-  a  a 

»■«  UJ  nJ 

ii  •  • 
Jx  X 
bJ  UJ 
0  0  3 

r*z  z 

iOhh 


*  <>- 

w  Q.  ^ 
*«  J3 
*3«t  * 

*“  I  H 

ft  to 

ft  ft  l~»  I 

*  r*  *■'« 


* 

* 

* 

«  -> 

*  -« 

*  ♦  „ 

*  X  *3  Z.O 

*  r;  *  *•  >-  uj 
**-*-»  ii  * 

^  ^ 

*  J  *  *  HI  „t 

*  <  *  w  w  _ 

*  2  #  t-i  z  < 

*“  *  >-*  *  I  V  x 

-  *  *  u.  *  a.  X 

x  z  *  *  •  -O 

s^ss:*sjl  * 

*  J  #  M  *)  . 

*  3  ft  -  .  .  0T? 

*^JZp«>iUJUJ?” 
*<»||  .>-22  ^ 


< 

I 

3 

«h 

(V 

V 

ft 

3 


x 

<t 

tv 

ui 

• 

tv 

< 

I 

X 

-I 

c 

I 

3 

fO 

>- 

* 


X 
<t 
UI  ■ 

ft 


r~ 

a> 

<v 


H<5  HH 
II  -  ||  || 
^Zhd 
X 

^  u>in  o 


X) 


co  f*- 
o  o 
IO  *o 


3U.ll. 

3  H  H 

O 

«-l 

*o 


rj  •■<•**  o  ft  •u.z 
fOjiOiO#  *  ^ 


o  <  o  o 
3  0  0  0 


tv 

o 


»  * - -  3  Z  Z 

* . *£x *03 

*  *  *  >-  j  u 

* 

«o  in  4- 


o  r*> 
>-  »o 

o  o 
o  o 


-*  II 
to  f-t 


~  o 
3  X  »- 
*  X 
—  3  O 
*—  *-*  K  10 
C  >-  uj 
X  #  X  -% 
X  ***  o 
3  -» 

•  3  • 
t~*  o 

*—  •  uj 
X  3  • 
<  aJ 
lO 

UJ  _ 
*  ►-  Q 

-*  <*  z 


l»> 


< 
Z  I 


II  w 
3  Z 


•X 

UJ 


O  Z  <  Ii.  u. 

3  >-  I  _ 


O  ^ 


X> 


0  3  fO 


h> 

3 


O 

O 


o  o  o 


^HHsHsS??sf3o3|S|S|SS3a252Si2S°- 

A-48 


i^>wacM 


TR  5341 


#  * 


* 

« 

# 

* 

ft 

« 

* 

* 

* 

* 

* 

» 

* 

« 

» 

* 

« 


* 

* 

« 


* 

* 


# 

# 

*  h- 

* 

*  li_ 


*  O  ft 

*  ft 

*  z  ft 

*  o  « 

*  *H  * 

*  *-  * 

*  3  ft 

*  z  # 

*  3  # 

*  u.  ft 

* 

#  < 
ft 

ft  « 
ft 

ft  J) 
ft  ►-« 
ft 

*  Ul 

*  Ul 
ft  IT 


* 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 


ft 
ft 
ft 
ft 
ft 
ft 
ft 
ft 
ft 

#  x 

* 

*  ~  i 

*  3  *  -*  ♦ 

*  ►-  *  **  -» 

*  #h)“ 

*  *  II  II 

*  ft  -ft  M 

ft  ft  ft  h-  ►- 


o  3  u 


CM 


»0 


z  • 

I-  * 

<  Q. 


X  U. 
+  < 
CM 

*-  -I 
II  -J 
•n  < 

h-  3 


-»  Z 

a.  - 

•  w  -4 

t-  II 

Z  ►ft 

u. 

<  -ft 

-t 

-I  O 
-I 

<  o 

3  3 


-»  3 
3  - 
► 

IQ  ~ 

**  V 

>-  ft 
ft  ~ 

—  3 
3  « 


3  —  <f 

•  <  I 

CM  I  ~ 

- 3 

>~  3  - 
ft  »  *h 


■)aa 

4hh 

I  »  - 

3  »  ► 

•  t  t 

O  O  H  w  w 
Q  O  —  MM 

•  •  *-  r  x 

OOi<UJ32  » 

II  II  •  ~  II  II  «  • 

OHH  ||  ||  m 

3  ►  *Tm  >  »  m 

•io  ^ 

O  *  »*ftO  *  »aicm 

II  ***11**** 

'"'“''''O-^'-'^lOlO 


3. 

#  *  * 

o 

*  * 

« 

«  « 

»  « 

*  # 

«ft 

#  # 

H 

*  * 

* 

*  « 

*  » 

w 

*  # 

H 

ft  IO  ft 

X 

ft  -  ft 

CL 

*  CM  ft 

♦ 

*  •>  ft 

ft  -H  ft 

♦  -  ft 

ft  O  ft 

IO 

«  II  ft 

ft  3  ft 

ft  ft 

*  *  ft 

M 

*  —  ft 

X 

*  3  ft 

o 

#  *  ft 

* 

*  IO  ft 

ft  ~  ft 

ft  ►-*  ft 

*  X  ft 

M 

*  a  * 

«  « 

to 

*  -  # 

*  z  * 

M 

*  O  ft 

X 

ft  ►-•  ft 

CL 

ft  1-  ft 

+ 

ft  3  ft 

-  X 

•J>  13 


x  o  -  X  x  o  o 

<3  3  3  O  3  Q  3 


it 

<Q 


3 

o 


CM 


X 

CL 

ft 

r 


>-  3 

II  •  z 

-*  3-  3C 
-4  -O 
"  t  I- 

«r  —  u 

>-  »-ft  ct 


*  l:  * 
ft  3  ft 

*  X  ft  ~ 

*  ft  X 
ft  *-ft  ft  < 

*  X  ft  ~  ft 

*  X  ft  -1 

*  ft  »  -* 
ft  h-  ft  -I  - 

*  *“•  ft  >-»  ~4 
*3*1- 

*  -•  ft  <t  IQ 
ft  -I  ft  ft  I 

*  X  ft  -*« 

*  £  ft  H  *-> 
ft  ft  V 

*  *HO 

*  uJ  ft  ^  3 

*  K  ft  (M  O 
ft  <  ft  I  • 

*  J  ft  <  -4 


13 


*  3  ft  II 

*  U  ft  — 
ft  -J  ft  -H 
ft  <  ft  % 

*  3  ft  -ft 

*  ft  - 
ft  *  n 

*  *  I 

*  *  *  < 


ii  ii 

'►IOh 
-ft  m 
-  it 
ftoi n 
IO 

* 

X  o  o 
<  (D  3 


o  — 

Z  3  * 

•  I 

-*  O  «I 
II  II  II 

- - 

3  3 
t  »  » 

M  H 

o  - - 

X  /Q 
O  <  X 
3  Li  <t 


ft 


t  t 
X  X 

<t  < 

ft  * 

3  * 
-  X 
—•  <4 

~  u 
’  X  -» 
►  < 
Hal  I 
II  II  UJ 
SC-> 
3  X 
t  •  H 
If)  Hft 
O  —  3 
£  3 
3  <t  <t 
3  Ul  «3 


4 

<Q 

Ul 


< 

<M 

Ul 

% 

X 

« 

Ui 


ft-  « 


J 

o 

< 

X 


E 

ft 

J  I 
a  ii 

>4  CM 
3  X 


X 

< 

IO 

aJ 

■k 

I 

< 

CM 

UJ 

» 

CM 

X 


UJ 

Q 

<  X 
X  ft 

CM 
J  X 
3  II 

<  m 
3  X 


CM 

o 

-t 

in 

if) 

•O 

o 

o 

o 

3  3  3 


cM;M:M<M;S<MiM$°K>oio^oio^oio^O***!**!2;25$5SS2{,2;,'‘n'^,  ^  ° 

************  ***>**** 


A-49 


TR  5341 


«****• 


ftftftftb 

3 

3 

• 

» 

* 

3 

•ft 

3 

*  Q 

H 

* 

H 

* 

*  UJ 

w' 

H 

H 

#  Q 

lO 

*— » 

<n 

ft— » 

*  UJ 

X 

K) 

X 

IO 

*  UJ 

< 

I 

< 

X 

*  UJ  z 

♦ 

4 

+ 

4 

#  X 

o 

+ 

♦ 

*  H  UJ 

3 

ft** 

-ft  3  3  — 

o 

•  X 

• 

-ft  3 

3  -•  «  3 

o 

*  z  4 

«o 

3  - 

*•1-4  H  » 

• 

#  o 

^  M 

M  w  w  H 

4) 

#  cl  in 

1-4  *-* 

—  CM  CM  ft* 

\ 

*  -ft  3  Z 

3 

—  (M 

CM  I  X  CM 

•ft** 

*  K>  O 

• 

* 

CM  X 

X  4  <  X 

3 

#  O  *-< 

H 

H 

X  4 

4  «  »4 

* 

*  »Z  1- 

co 

w 

<r  * 

*  O  O  » 

1-4 

«  jJ  Ul  4 

z 

(M 

*  o 

o  n  o  o 

** 

*  J  CL  3 

z 

» 

I  O 

O  O  O 

D  in  -»  in  d 

CM 

*  •  UJ  <  3 

1-4 

<!> 

4  O 

Q  Q  If)  — 

IT  •  3  •  • 

X 

*  0X0 

* 

*  • 

•  •  •  3 

•  -4  ft  O  IO 

< 

•  X  CL  3 

Q 

o  io 

K)lO  O  * 

H  I  H  +  | 

* 

*  UJ  O  3  < 

3 

H 

o  \ 

\  +  1  M 

♦  -ft  w  -«  — 

o 

*0  1-41) 

X 

•b 

•  •*•»* 

^  w 

^  3  CM  3  3 

o 

*  L 

O 

H 

H  3 

333CM3«X-- 

• 

*  3  3  ft  • 

1-  CM 

% 

H  • 

-  *  *  I 

ft  H  4  H  H 

H 

X 

«  z  ui  in  < 

in  cd 

H 

i  •-• 

o 

M  H  H  <  M  ^  *f  ww 

•H 

UJ 

*  »—»—»— 

•ft# 

Q 

i 

w  X  -ft  X  I 

♦ 

Q 

#  CL  <  Z  UJ 

3  Z 

CM 

3  m 

•!\II  I- 

14  344 

Z 

*  UJ  3  iU  13 

x  ►-* 

■« 

-  X 

>014  <T<  *  *  *  * 

3 

H 

*  H  3  — 

4 

% 

*1  4 

\<l  *  *  * 

I  OMOD 

* 

* 

*  in  i.  o  in 

* 

—  1 

i  oob-ioO'-'Qn 

H 

♦  -1CHU 

in  in 

*■# 

I  — 

3 

'*3  0" 

3  •  I  •  • 

>* 

*  in  o  u.  *- 

3  D 

* 

<  3 

* 

3  -  -X 

•  ^  <  cm  in 

X 

* 

«  O  U.  U.  4 

O  O 

H 

#  - 

M 

•■  ID  C\J  < 

j-  ♦  i  i 

4 

2 

#  O  3  3 

X  X 

* 

O  *- 

*-l  1  1  ♦ 

•f  •*  ^  * 

* 

* 

#  X  3  O  3 

I-  I- 

H 

Q  *- 

CM 

W  «ft-% 

•*  t  n 

o 

>- 

*  1-  X  O  O 

3  3 

%-*  o 

•  I 

X  I  3  3  3  3  *  -  «  - 

o 

* 

#34  3 

X  X 

h  r- 

CM  4 

•  H  H  M  H 

• 

X 

«  y  o  <* 

-5  iO 

♦  ♦ 

♦ 

|  H  H  H 

H  w  w  w 

CM 

* 

*  3  •-•  U 

1—  h 

«k 

<»"*  ^b 

ft  *— ft"  '  «— * 

—  *-  i-  «  -ft 

♦ 

Q_ 

*  CL  1-  V- 

H  H 

X  o 

3  3 

33I-»-I-Kmhmm333 

** 

UJ 

*  ui  z  mo 

<J  o 

<  h 

*  •» 

•ft 

*mhmh2ZZZ » » » 

3 

t- 

#  hlJHO 

H  H 

|Q  >-«•—« 

UJ  O  — 

ft  O  K  t— 

K)  M  M 

I  ~  z  z 
*  o  z  z  3  o 

<  -  -  I  I 

-  •  *-!  *4  ||  || 

3  0  II  II - * 

QIiJhTTT 

•  »  1 « 
CL  4)  Ahm 

k  in  j)  •  » 

0  —  - - 

4  U.  O  O  *4  *-* 
U  «-*  O  Q  * 


y-  Z  Z  TTZOODO-**-*!-*  - 

-'—0333#*#*  •  -  -  n 

»-►-#  *  #  tOOOONfli- 
►-••-•OOOOOOQQft'—ft-*- 
Z  X  3  Q  3  3  •  •  • 

03  •  *  •  •tomioiojeji.tz 

I  llOlOKIlOt  I  I  I  I  I  o 

II  II  II  II  II  II  II  II  II  II  II  II  II  II 


JE  _«  Z 

I  I  3 

II  II  II 

-ft  —  -»  o 

lb  »  ft 


in  *  in 

*—#!—• 
in  #  3 1 
X  #  3  ' 

J  t  i  I 


- - ' - -  —  — - - -  —  Q 

**ft-tiM<\IC\IC\lfOrO’OfOSha'4’3-  O  Z 

*  *  *  *  *  *  *  ><<X<XXOiJ 


I^OUlH 

3 

O  >-  U-  lL 
X  3  O  O 
i  3 
»-  o  in 
< 


■*##### 


uuuuuuuuuu 


<o  m>  «o  n  <o  3  bOvOiOr-r^r-h-r-i^i^r-r-r^co*^) 


^cu'o^-bn^r^^a'o^ 


THIS  bUbrxOUTlNE  IS  CALLED  BY  START  OR  DIFEG 


COmMON  AU0#20)  *  ALPHA  (4)  ,T<20) 

DIMENSION  Y<4,20) ,YN<20) »B<4) ,81(3) »B2(4) ,FN(20> 
DOUBLE  PRECISION  A , ALPHA »b*Bl»B2»FN,H»T»Y»YN 
DATA  1 BETA/0/ 


TR  5341 


■> 


o 

Q 

o 


CO 


o  o 
Q  Q 


o  o  o 
,-J  •  • 
O  O'  CO 

•  +  \ 

o 

CO  o  o 

Q 

X  D  Q 

in  o 

o  o  o 

r-  q 

D  •  • 

•  o 

O  rO  O' 

o  •  o  •  \  + 

+  IO  Q  K)  ^  O 

osm  +  ioo 

C  ^  (\i  o  — ■  o 

0*0  •  Q  «t  • 

•“NOIIO 


2.  o 

•  Q 

H  •<  H  • 

a.  ii  ot-  -h 
UJ*-»  II  UJ  II 
I-  •  *  2)  O 

l/l  *4  *-•»-•  |-  -4 

II  o 
*  a 


o  o 

K)  <  +  •  0.  v. 

O  o 

\ioioj« 

o  o 

o.  a  v  <  ro 

•  • 

IO  J  o  ~  *  ~ 

to  m 

<T  »K)  O  «X 

\\o 

«  *  *o  ~  a  x 

o  o  a 

x  o  \  <  o  a. 

a  o  o 

a  o  ~  x  •  _J 

•  o  • 

JOIOUt  < 

~4  •  fO 

<  •  —  _!+  + 

♦  ^  V 

+  d-  <  «  o  o 

0+0 

o  <  x  v o  o 

o 

a  o  a 

Q  o  a  o  o  o 

-4  OL 

o  a  o 

o  a  _i  a  -  •' 

•  O  • 

•  O  <  CD  +•  J- 

*- 

CM  •  *4 

CM  *  +  •  CM  CM 

◦  in 

HIO  + 

~t  *)  O  JJ  V  V. 

*-  X 

WO 

WQ\''-' 

— - a 

-'^o-dcM 

o 

CM  CM  CM  • 

CM  CM  »CM  — '  *<*■ 

CD  o 

«w(\| 

"  ~  CM  ^  4  < 

o 

<  O  <  <  *+ 

<  <  H<  x  X 

-  *o 

i  3  x  r  Nk 

I  X  X  I  X  X 

o  - 

a.  o  a.  a.  ~ 

X  X  X  _*  _J 

o 

— 1  •  _J  _l  <M 

J  _J  cm  .  J  <r  <t 

•  o 

<  CM  <  -X 

<<-’<## 

I-  CM 

#  +  *  #  < 

*  #  <t  #  o  o 

CD  *  O 

o  o  o  x 

o  o  i  o  a  q 

•  o  o  a 

-4  J  H  Q  □  1 

x  3  x  a  o  o 

OQlf) 

-oin-oo  jo 

o  o  J  o  «  • 

o 


•  <MO  »H  •  •  <  O 

OIQ  ^02)2101  I  cf 
II  II  II  II  II  II  II 
O  ^  ^  ^  ^ 


If)  CM 
II  II 


a 

«  •  O'  in  _l 

I  lOH  |  4 

II  II  II  II  II 


-Tu.OH(gniOHHN(\JNO 
JJ2  mSDOSlS 


»4  -»  -i  <M  CM  CM  CM 
O  -O  D  X)  □  lD  O 


o 

o 

O  ul 

II 

<  o 

*-  I— 
UJ 

D  O 

*-<  o 


•-i  o 


C*  I—  K) 
UJ  o  * 
•  »H 
II 

a.  x  h 

HUlaJ 
II  I—  .J  If) 

I- 

Ul  —  ~ 

)u.u.o 
<-**-•  -*  D 


u. 

* 


o 

o 

CM 


o 

o 

fO 


O  O  If) 

O  -4  rl 

*  *  J- 


iM,o^’<nor»0O'o<-4CM'O4'if)<ON‘CO9i>o>^^'O4-in«or>2>O'O'4iMio4’tf><or*> 
•4«4r4.4.4.4~4.4CMv\|CMCM>'MtMCM.MCMCM«)i0<0''>«0i0  'ifO/OfO****  *  zT  *  3 


A-51 


ft 


420  DO  450  I=i*KSTEP 


TR  5341 


ft  « 

ft  ft 

* 

ft 

# 

ft 

* 

ft 

ft 

ft 

ft 

ft 

ft 

ft  -V 

ft 

ft  o 

ft 

ft  IM 

* 

ft  — 

» 

ft  3 

*  2 

ft  o 

ft  3 

ft  J 

ft  ►* 

ft  - 

ft  1- 

ft  ~ 

*  < 

ft  © 

*  X 

ft  CM 

*  H« 

ft  - 

ft  X 

ft  o 

*  o 

ft  CM 

«  tr 

Jl  — 

ft  a. 

ft  u 

«  a. 

ft-  - 

«■  < 

ft 

* 

.*  o 

.*  UJ 

.ft  Of 

ft  o 

ft-  * 

ft  < 

ft’O 

«  a. 

ft  CM 

« 

ft  — 

*  >- 

*  ao 

ft  CO 

«■  * 

ft 

.*  ~  X 

ft  3 

ft  o  K 

<■*.*« 

.ft  CM  O 

Z  ft  •-« 

ft  »z 

•  ft  I— 

ft  O  X 

3  ft  2  UJ  ft  CM  - 
O  ft  UJ  2  ft  —  UJ 

U  ft  2  m  ftX  > 

in 

M 

»  ft  O  1-  ft  *  < 

<c 

O  ft  X  3  ft  -»  X 

4 

—  3 

-  ft  X  O  ft  o  - 

H 

3  — 

ID  ft  UJ  X  ft  CM  X 

ftu 

o 

-  >- 

*  ft  CO  ft  -  » 

3 

b- 

2  * 

a  *  x  3  «  o  « 

U.  ~ 

*  #  -H  in  ft  n  i- 

< 

O 

ft  3 

x  #  a  *  —  ui 

(J>  -> 

•*-* 

•  ft  1-  U)  ft  <  CD 

in 

X 

2  «t 

<t  ft  <  X  ft  s 

3 

% 

—  XX 

-ftXOftZZO 

o  o  o 

< 

X  2 

CD  UI  X 

ui  *  z  «  o  o  a 

2  3  3  3 

3 

UJ  X 

*  1-  -1 

o  *  in  *  hmio 

»  •  •  • 

♦ 

t-  * 

x  in  < 

<#w>-*inino 

HOQOU 

in  z 

Z  +  2  X  1  J- 

IfthUftWHlO 

II  II  II  II  3 

H 

x  * 

ft  <  ft  O  O  • 

-3  ~  ~  ~  £  o 

2  2  — 

Ui  <3 
3  UJ 

2.  • 

*-  < 
Z  « 

o  x 

U  H 


—  II 
Z  'I 
u. 

U.  o 
>0 
J  4 
3 

<  o 
o  3 


H  II  M  p 

•f  »-i  —  CM 
<M  ah 


II  II  — 

HlZ  UJ 

>-  3  _ 

OOlIZ^iilOO 
NN-»M(£KHllO 

^  a  h  h  3  >/)in  h- 
—  2  H-  *  *-<  O 

OOZOUIIO-O? 
03>-Uir^Qlll9il 


in 


ui*30ftuiuiN, 
2*OOftir3:< 
~i  ft  O  J  ft  X  X  I- 
t—  *  _l  _l  *  UJ 


2  ~i  3  -«  a 


II  o  —  —  — 


2  O  UJ  II 

o  II  3  H 


-t  o 
II  o 
3  II 


O 

o  in 

O 

o 

<0  vO 

r* 

o 

4- 

ah  4 

4 

n 

3  ft 

<<*UUS£hD 

3  o  a  O  ~  z  ~ 

O  ft 

UU*  J  J  II 

M  H  f»  0  H 

X  ft 

ft  33  3  <  UI  CM 

—  ft-  **  H  — 

3  ft 

ft  3  3  1-  > 

3  2  3 

3  ft 

ft  O  O  <  <  O 

O  O  O  O  O 

in  « 

**#30313 

3  0  3  3  3 

-I  CM 

(JUJU 


ioo>oH(\in^inDMo»uH!M 

4-^inin<nimnintnininu)va>o<o 


^(\iiod-nnf^»J*o^cM^^'n«r^»j'o 


A-52 


TR  5341 


<n 

r( 

o 


o 

3 


J 

o 

o 


*  « 

« 

« 

# 

* 

* 

* 

# 

* 

* 

#  z 

* 

*  o 

* 

*  *-• 

« 

*  H 

« 

*  U 

* 

#  3 

* 

ft  3 

* 

«  jJ 

« 

#  X 

* 

* 

* 

*  1/5 

« 

#  3 

« 

#  oJ 

* 

«  xl 

* 

*  Z 

* 

« 

.* 

*  — 

* 

*  < 

• 

#  — 

* 

*  a 

* 

.*  X  Jk 

♦  UJ 

* 

iZ(5  h 


xl  xl 
3  3 
Z  Z 
— *  <-« 
H  H* 
Z  Z 
o  o 
o  u 

vO  f*» 


(3 

-J  "* 

O  II  X 

O  — <  X 

II  o 

S  flZ 

I  HX 

o  — 1 

Z  o  x 
X  3 


O  .J 
3  3 
II  2 
X  >-i 

or  *- 
O  2 
Z  O 
X  3 


*  * 

*  oc  # 

»  aJ  # 

*  X  * 
#»-  * 

*  UJ.* 

*  X  * 
.*  m  * 

*  # 

*  uj  * 

*  UJ  * 
*l/»  * 

*  *• 

*  O  * 

.* *-  .* 

•  * 

*  3  * 

*  uj  * 

*  i/i  * 

*  3  * 

*  # 

*  l/>  * 

*  * 

#  2E  * 

#  X  * 

#  o  * 

#  Z  # 

#  * 

*  z  * 

»  X  * 

#  3  # 

*  3  # 

#  O  * 

*  3  # 

#  * 
«  « 
*  « 
*  *  * 


o 

CM 

w> 

o 

f\l 


o 

3 


I 

X 

* 

X 

cr 

o 


X 
o  -» 
II  u. 
x  H 


*  *  * 

*  # 

♦  * 

*  * 

*  # 

#  * 

*  4 

*  * 

*  —  * 

*  <t  * 

*  *  * 

*  *"»  * 

*  «  * 

*  ►-  * 

*  UJ  # 

*  rr  * 

*  l  * 

*  in  * 

*  •  * 

#  ~  # 

.*  +  .# 

*  -<  * 

#  —  *■ 
*  *  * 

«  *->  * 

»  I  » 
«  «"  * 

4  *  * 

*  *  * 
*  ~  * 
*-t\l  * 
.*.*  * 
*  >  * 
*  <  * 
*  *  * 
*  <  * 

*  P-  # 

*  UJ  * 

*  tn  ♦ 

*  ♦  # 
*  <  <* 
*  *  * 
*  ~  * 
#  <  * 

*  P-  * 

#  UJ  * 

*  CD  * 

#  +  # 

*  \n  jt 

#  •  # 
*  «-  * 

*  I  * 

*  *•«  * 
*  -  * 

*  ii  * 

*  ~  * 
*  «  * 
*  -  * 

*  o.  # 

*  X  # 

*  UJ  * 

«  * 

*  * 

*  # 

#  #  # 


X 

* 


< 

« 


1  % 

*  1 

<  Q 

4  m 

• 

*  T 

-1  X 

~  .* 

<  — 

♦  3 


z  (-1  0.  ° 

*  W  *  3 

H  Ql  al<  ® 

II  "  5  K,  J 

z  *  —  z  uj  7 

«.  th  a  + 

3  •  p-  -»  UJ  — 

II  MZIIDH 

~i  O  o  »■*  2  * 

*n  ‘ 


X 

* 


< 

4 


xl 

X 

I 

o 

Z  3  O 
*  X)  3 

H  *  • 

O  UJ 

II  II  3 


o 

Q 

© 


II 

3 


II  O 
m  3 

>0 


X  Z  O 

jj 

hZU  ^ 

wOll  3  2  110—^3 


J  UJ 
II  3 
-  2 


II  II 

r-«  ~i 
CM  -O 


o 

•O 


o 

D 

K) 


«.  H  J  O' 
H  Z  J 

~  o  <  o 

u  UU  3 


in 


M  z 

w  O 
3  3 

O' 


o  o 

3  3 


3  3  3 


3  3  3 


A-53 


UO  1H  K=1»N 

p ( 1 »  J ) sP 1 1 » J ) 1 1 »  K ) ( K » J  i 
m  CONTINUE 
13  CONTi-'lUt 


.<OKM(AH)  .GT.  (.1)  »  EXP(A)  =EXP  t  A/2*£m )  ( 2**M ) 


TR  5341 


« 

« 

* 

* 

« 

* 

* 

* 

* 

« 

* 

* 


►-  X  * 

Z  -»  ft 
O  UL  ft 

3  *-•  * 
CM 


*  « 

* 

* 

* 

* 

* 

* 

# 

* 

* 

* 

« 

» 

# 

* 

* 

« 

* 

* 

* 

« 

* 

* 

« 

.* 

#• 

* 

« 

* 

* 

♦ 

.* 

* 

# 

« 

# 

* 

* 

«■ 

* 

.* 

* 

« 

* 

* 

* 

« 

* 

# 

« 

» 
ft 

ft  2 
ft  - 
ft 

ft  II 
ft  *-• 
ft 

ft  4- 

•  rvi 

ft 

ft  o 

ft  ft  3 


z  o 

-  3 
H  • 
II  O 
3  II 

in  3 1 
CM  *  l 

t-« 

O  « 

a  tn 


Ui 

3 

Z 


CL 

ft 


QL 

♦ 


at  z  z  z  * 


z  z 


Z  UI  II 
o  3  * 

a  z 

h  m 
z 

o  o 

3  3 

«n  ah 

CM  cm 


ii  ii  ii 

*-«  3  _) 

-o  o 

CM  CM  cm 

O  O  O 
3  3  3 


CQ  ui 
II  3 


3  jJ 
3  3 
Z  Z 


II  II 
>-<  o 


o 

'  o 

• 

o 

ii 


z  z  o 

»  •  Q 


o 

3 


3_hhh!M33 


-•  Z 
—  O 

n  o 


z  z 
o  o 

3  3 


<M 
K>  >0 

O  O 
3  3 


0.  V 


UJ  3  Ui 
3  3  3 
z  z  r  ui  z 
«  >  <r  <m 

hl-HOS 

z  z  z  r  »-  r 
5  o  o  ii  ui  ii 
3  3  3 1  ae  r 


ii  ii 

H  3 

i  -t  iO 

m  m 
o  o 

3  3 


0  3  3  0 

II  3  3  K) 

—  Z  Z  3  Z 
3  <m  _  O  >  (K 

-•  Z  Z  Ih 
w  O  O  O  II  3 
X  3  3  3  r  X 


O'  35  r- 
CM  CM  CM 


CM  3 
<0  -O  (O 


O 

CM 


in  4- 
in  n 


o 

4 


3  3  3 


jiOH«ifl#iniONSJ,oHaioafli)M)jioHiMC)ai)OM)>oHiM<)tin 


A- 54 


TR  5341 


«  • 

* 

# 

# 

* 

# 


*  *  * 


# 

* 

* 

* 

* 

*  UJ 

*  O' 
ft  J 

*  X 

* 

*  a 

*  UJ 
X  ft  ft- 
•  *  2 
~  *  uj 
o  * 

<U  #  W 

~  *  r 

K  ft  C 

•>  ft  ti 

~  *  K 

o  *  < 

<vj  *  _j 

~-  *  3 

tD  ft  u 
-  ft  _J 
«’"»*< 
o  #  o 
w  *  CM 

«  ft  _J  • 

<  -  ft  «  £ 

*  >-  #  + 

»-  •  #  .  z* 

*  ~  ft  —  w 

3  O  ft  v.  + 

*  «t  *  •  >- 

>-•*►-  o 

►  O  ft  w  o 

X  <M  ft  O  ^ 

X  «  *  UJ  U)  l| 

*  ^  I  ^  H 

°  *  X  X 
2  ►■»  *  in  v 
U-  00  ft  aJ  >-  ~j 
‘3  «  *  »-  a] 

O  ft  < 

J  J  *  J  a 

—  #  O  UJ 

*■•  0-  *  O  3  • 

*“  *  3  3  iD 

—  J  #  <r  <1  » 

^  J  *  y  u  u 

X  o  # 

0  3  ft 

-3  3  ft 

«/)  J  *  »  «  § 


# 
# 
ft 
ft 
ft 
ft 

—  ft 
3  # 
*-  * 
t-  # 
ft  # 

—  ft 
“3  ft 

—  ft 
t-  • 

♦  * 

•  * 

•ft  * 

II  * 

—  « 

•ft  « 

“  # 

<9  ft 

ft 

* 

uJ  ft 

2  * 

►ft  # 

U.  ft 

3  *  ~l 

■3  *  r 

*  3 

ft  I- 

ft  ui 
*  ft 


* 
ft 
ft 
ft 
* 
* 
* 
* 
» 
* 
* 
ft 
* 
* 
ft 
# 
* 

* 

ft 

ft 

« 

ft 

* 

# 

* 

ft 

ft 

» 

# 

ft 

ft 


■3 
3  X 
X  J 


<-><->0  0  0  0 


"tc\i>Qd.tn,0l^xl3,o 


*  ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 


ft  ft 


ft  ft 
ft 
ft 
ft 
ft 
ft 


ft 

*  UJ 

*  cr 

*  jj 

*  l 
ft 

*  a 

*  UJ 

ft  b- 

*  2 
ft  UJ 
ft 

»  W 
*  2 

*  O 
ft  i-» 
ft  K 

*  < 
ft  -J 

#  3 

#  O 
ft  _J 

#  < 

^  *  O 

~  ‘  ♦  J 

5?  *-  ft  -I 

*M  -  ft  «* 

W  T*  M 

*“  U.  ft  • 

*  *  ft  ^ 

^  o  r  *  . 

;-iva»K 

■2  <  *  3  ft  w 

U.  I  ;*  4  «  u. 

»  a.  —  »  * 

2  3  >-  X  ft  uj 

»  <  »  #  x 

>“  »  ~  X  ft  K 
X  -*  O  o  ft 

2  O  r\l  H  ft  tn 

u.  cm  —  in  #  uj 

U-*2>-**i- 

o  u.  U  ft  < 

UJ  CM  UJ  ft  3 

z  ~  z  <r  *  z> 

“•  <  o  Q.  ft  (J 

•”  ft  3 

3  Zinm  *  < 

o  O  2  J  #  o 

CK  *  UJ  CD  ft 

JilO# 

-3  O  i-i  o  ft 

^  U  O  J  ft  ft 


(M  + 

ft 

ft 

in 
x  ft  « 


—  * 
K  ft 
ft  ft 

'*»■  ft 

►1  ft 
—  ft 
I-  ft 
ft  ft 
•  ft 
•ft  ft 
ft 


ft 
ft 
ft 
»  ft 
ft 


UJ  ft  ft 

2  V  . 
►ft  o  o 
K  O  O 

3  >ftH 
O  I  I 
(T  II  || 
a  i — > 

3  a  *ft 

in  s. 

>-  2 
>:  a  u. 

CO 


ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 


a  uj 

UJ  2 

-J  »H 
JOlL 
<  •  UJ 

u  UJ  3 


ft 

ft 

ft 

ft 

*  2 
♦  X 


*  3 
*»-3 
ft  aJ  2 

*  ft  ft  ft  x  a 


u  u  u  o  u  o 


H  ^  lO  J  JI  |s  B  JIO 


CM 


• 

ft 

ft 

ft 

ft 


ft 

ft 

ft 

ft  uj 

ft 

«  a: 

ft 

ft  UJ 

ft 

#  X 

ft 

ft 

ft 

ft  in 

ft 

*  2 

ft 

ft  o 

ft 

ft  >~« 

ft 

ft  H 

ft 

#  < 

ft 

ft  h- 

ft 

*  3 

ft 

*  a 

ft 

#  5l 

ft 

#  O 

ft 

*  U 

ft  * 
ft 
ft 
ft 
ft 
ft 
ft 
ft 
ft 
ft 
ft 
ft 
ft 
ft 
ft 
ft 
ft 
ft 


ft 

ft 

ft 

ft 


Uj 


I—  o 
-  CM 


•  <r 
< 


*  3 

*  3 

*  « 

# 

ft  cr 

*  uj 

*  t- 

* 

* 

« 

*  • 

*  l- 
« 

*  u. 

*  o 

b-  ft 
-  *  2 

~  *  O 

O  ft  -I 

ft  I- 
#Uf- 

*  ■-  o 
*30 

*  U  -4 

2  • 

*  <  II 

ft  -V 

*  m 


CM 


>-  o  „  _ 
z  -*  ft  *-• 
u.  <n  * 

4  »-« 

u 
UJ  3 

2  ir 

*-«  a. 


t 

*  4 
ft 

*  X 
ft  -1 

*  X  • 

-  .  *  K  3 

3  3  ft  <  • 

‘  *  *  U 

ft 
ft 
ft 

*  ft  ft 


ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

I-  ft 

*  ft 
•  ft 

O  ft 
o  ft 
•ft  ft 

I  ft 

II  ft 
-*  ft 
•ft  ft 

•  ft 
ft 
ft 


3  J 
a  jo 

3  3 
3  O 
on  -j 


<  * 

* 

« 

UJ  ft 

2  ft 
►•*  ft 
U.  ft 
UJ  ft  2 

3  ft  X 

*  3 

*  *-3 
*32 

*  #  r  a 


<->  3  U  3  U 


^  cO  *  />  ,0  f».  s 


'•0  » 


A-55 


SUbrtOUTlNE  PrfINT(H»T#Y,N#I> 
OOUfaLE  PRECISION  H,T»Y(20) 


TR  5341 


*  *  * 

* 

#••••# 

* 

• 

•  * 

* 

« 

•  » 

* 

• 

•  * 

* 

• 

♦  * 

• 

• 

*  • 

* 

» 

*  * 

*  b- 

• 

*  * 

* 

• 

*  Q  * 

#  >- 

• 

*  UJ  * 

* 

• 

•  X  * 

# 

• 

*  in  ijj  • 

#  cn 

# 

*  s:  i-  • 

*  t~ 

• 

*  _i  _j  • 

*  J 

ft 

*  2  <  * 

*  D 

• 

*  2  * 

#  1/1 

* 

#  X  O  « 

*  i 

* 

*  O  * 

#  X 

* 

*  in  • 

* 

* 

*  UJ  2  • 

«  u 

* 

*  .3  *4  • 

*  ►— 

• 

•  4  4  « 

*  4  l/) 

• 

*  X  X  * 

#  *“<  K 

• 

«  UJ  X  * 

•  04 

• 

*  >-  X  UJ  * 

*  U»  X 

• 

•  cd  m  * 

*  i  x 

• 

*  fl  73  * 

#  tr  o 

• 

*02  * 

*  U  U. 

* 

*  UJ  4  UJ  • 

*  t- 

#  ^ 

*  _J  X  * 

*  /  h 

# 

#  _i  in  h  * 

*  •-•  3 

•  X 

*4  1-  « 

*  o 

*  cu 

*02  >-  * 

*  h-  »- 

*  - 

•  UJ  CD  * 

*  3  z 

»  CD 

*  ••  ar  • 

*  O  i-t 

*  • 

*  UJ  UJ  UJ  * 

*  H-  x 

•  in 

*  2  -I  _J  * 

*  2  X 

#  H  ^ 

*  t-t  UJ  CO  * 

*  H 

*i)2 

*  J-  ~  4  # 

*  X  z 

*  * 

*  3  _J  — «  UJ  * 

*  Q.  * 

*  in  «-• 

*04  1  O  * 

*  O 

*  ■>  ii 

in 

#  X  2  —  4  * 

*  3 

*  x  -j 

2 

*  CO  *-<  *  _J  * 

♦  l-m 

*  &  * 

4 

*  3  <3  •  CL  * 

*  >-* 

*  -  -« 

«k 

*  in  »-«  4  uj  * 

•  XX 

* 

Z 

*  X  X  • 

*  UJ 

*  X  — 

*  2  O  UJ  • 

*  J)  oJ 

*  m  >- 

4 

*  o  i  m  * 

*02 

*  -■  — 

t  W  U  t-  •-»  # 

*  *-• 

*  X>  - 

»- 

*  in  r  * 

•  XU. 

*  •  y~ 

X 

*  x  »-  m  x  * 

«  O  ui 

*  in  « 

oJ 

*  UJ  2  4* 

*  U.  O 

*  -H  X 

> 

•  >  in  *-•  x  * 

«  oi 

3 

*  2  2  4  ID  * 

#  ^ 

M 

#HMKO# 

*  UJ  < 

*  CM  (\l 

*  4  2  X  * 

*  Z  O 

*  *  H 

UJ 

*  X  *-  o  x  * 

*  3 

#  K  • 

z 

*►-<2  0  * 

♦  HU 

*  -*  *  UJ 

►-« 

*  x  o  in  * 

*  21  W 

* - o 

*  *-  o  in  * 

*  -U  l/> 

*  V-  z  z 

o 

*4  2  X  *  2 

*03 

*  <tniHir 

o 

*  3l  4  4  h-  •  X 

* 

*  X  K  b-  3 

X 

*  *3 

• 

*  ac  —  z*-  o 

0) 

•  *  i— 

* 

*  3  x  o  ui  2 

D 

•  •  UJ 

*  *  * 

*  U-  *  O  X  U 

in 

<M  ^  or 

CO  CO 


ouou 


O  O  U  O  U  O 


•HlU’O-i’lOsOf^-lOCT1 


A-56 


1..  - 


INITIAL  DISTRIBUTION  LIST 


Addressee  No.  of  Copies 

ONR,  Code  427,  480,  481,  483,  485  5 
CNO ,  Code  OP-2  3T  1 
NRL  1 
NAVELECSYSCOMQY,  Code  03,  051  2 
NAVSEASYSCOMHQ ,  SEA-0 9G3  4 
NAVAIRDEVCEN  1 
NAVWPNSCEN  1 
DTNSRDC  1 
NAVCOASTSYSLAB  1 
CIVENGRLAB  1 
NELC  1 
NUC  1 
NAVSEC,  SEC-6034  1 
NAVPGSCOL  1 
APL/UW,  Seattle  1 
ARL/Penn  State  1 
DDC,  Alexandria  12 
Polytechnic  Institute  of  New  York, 

Professor  Stanley  Preiser  24 
Westinghouse  Electric  Corp .  ,  Mr.  J.  T.  Malone  4 
Consolidated  Edison  Company  of  New  York,  Inc. , 

Mr.  P.  Hsiang  1 
Mr.  R.  Casal,  Brooklyn,  New  York  1 
RDP  Inc.,  Mr.  C.  A.  Steele,  Jr.,  1 
Sperry  Systems  Management,  Mr.  James  Leong  1 
Babcock  6  Wilcox,  Mrs.  Anne  M.  Schwartz  1 
University  of  Illinois,  Urbana, 

Dr.  S.  W.  Joshi  1 
Clarion  State  College,  Dr.  Stephen  Gendler  1 


